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ABSTBACT 


This  thesis  presents  a  Nuclear  Targeting  and  Effects 
prograa  that  is  intended  for  inclusion  in  the  simulation  of 
Tactical  Alternative  Besponses  (STAB)  combat  model.  It  is 
presented  as  a  stand  alone  program,  written  in  SIMSCRIPT, 
which  can  be  easily  modified  as  a  subroutine  for  any  high 
resolution  combat  model  requiring  tactical  nuclear  effects 
simulation.  When  presented  with  a  group  of  targets  which 
are  deemed  suitable  for  attack  by  tactical  nuclear  weapons, 
the  program  will  select  units  to  fire,  select  proper  yields 
for  multiple  yield  weapon  systems  and  assess  casualties 
among  Armcr  and  Infantry  within  the  target  area. 
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I.  IHTBODOCTIOM 


This  thesis  presents  a  nuclear  targeting  and  effects 
program  that  is  intended  for  inclusion  in  the  Simulation  of 
Tactical  Alternative  Responses  (STAR)  combat  model.  It  is 
presented  as  a  stand  alone  program,  written  in  SIMSCRIFT, 
which  can  be  easily  modified  as  a  subroutine  for  any  high 
resolution  combat  model  requiring  tactical  nuclear  effects 
simulation. 

When  presented  with  a  group  of  targets  which  are  deemed 
suitable  fcr  attack  by  tactical  nuclear  weapons,  the  program 
will  select  units  to  fire,  select  proper  yields  for  multiple 
yield  weapon  systems  and  assess  casualties  among  Armor  and 
Infantry  within  the  target  area. 

1.  CAPABILITIES 

The  model  in  this  thesis  analyzes  the  given  targets  for 
destruction  by  nuclear  weapons.  In  so  doing  it; 

1.  Can  handle  any  number  of  targets. 

2.  Can  handle  any  number  of  nuclear  capable  firing 
batteries  each  of  which  may  have  any  number  of 
different  yield  rounds  available. 

3.  Can  keep  track  of  any  number  of  friendly  maneuver 
units. 

4.  Erovides  for  friendly  troop  safety. 

5.  Represents  targets  as  composed  only  of  personnel  and 
tanks. 

6.  Assesses  damage  to  each  individual  target  element  in  a 
stochastic  manner  rather  than  as  an  expected  value. 

7.  Assesses  damage  caused  by  blast,  heat  and  radiation. 


11 


B.  LIMIIITIOMS 


Seme  of  the  items  listed  as  limitations  are  such  only 
because  net  all  weapons  are  included  in  the  data  base  or  nor 
all  possible  target  elements  are  included  in  rhe  targets. 
Host  of  the  following  can  be  input  into  the  model  with  a 
moderate  amount  of  work.  The  model  does  not; 

1.  Dse  free  flight  rockets,  missiles  or  Air  Force 
delivered  munitions. 

2.  Use  strategic  weapons. 

3.  Consider  collateral  damage  avoidance. 

4.  Model  any  other  effects  such  as  NBC  warning  and 
reporting,  optimum  time  of  exit  for  units  caught  in 
fallout  areas,  crossing  fallout  areas,  decontamination 
or  medical  evacuation  and  treatment  of  casualties. 


C.  THISIS  EHETIEW 

Chapter  II  will  briefly  review  the  effects  of  nuclear 
weapons.  The  damage  mechanisms  will  be  introduced  while 
formulas  and  methodology  for  building  the  assessment  phase 
of  the  model  will  be  developed. 

Chapter  III  will  discuss  some  of  the  more  relevant 
definitions  the  reader  must  become  familiar  with  and  then 
will  shew  how  the  analyst  in  the  field  solves  the  target 
analysis  problem.  The  model  will  mimic  this  manual  solution 
in  order  to  arrive  at  the  same  solution  a  field  anaylst 
would  get. 

Chapter  IV  will  discuss  the  mechanics  of  how  the  model 
works  without  being  code  specific.  Some  simplifying 
assumptions  will  be  made  and  approximations  will  be  offered 
to  reduce  the  processing  time. 


f  .  -  ,  '  .  »  »  _  »  ^  1  .  .  .  .  A.  . 


Chapter  V  will  explain  the  code  of  the  model.  Each 
routine  and  function  will  be  discussed  separately.  All 
global  variables,  permanent  and  temporary  entities  along 
with  their  attributes,  set  membership  and  ownership  are 
explained  and  must  be  fully  understood  to  comprehend  the 
model. 

Chapter  VI  will  identify  possible  areas  for  further 
research.  Areas  include  a  method  for  optimum  assignment  of 
batteries/yields  to  targets,  standard  field  operations  taken 


during  and  after  a  nuclear  burst  and  inclusion  of  ether 
elements  inxc  the  target  area. 


II.  EFFECTS  OF  80CIBAR  BEAPOHS 


A.  I1TBCDOCTIOH 

A  nuclear  explosion,  like  any  conventional  explosion, 
results  frca  the  very  rapid  release  of  a  large  amount  of 
energy  in  a  small  space  -  an  energy  density.  In  simple 
mechanics,  the  liberation  of  energy  is  manifested  as 
pressure,  heat  and,  in  the  case  of  nuclear  weapons,  ionizing 
radiation.  The  damage  mechanisms  of  nuclear  weapons  are 
refared  to  as  blast,  thermal  radiation  and  nuclear 
radiation.  Damage  caused  by  blast  is  further  divided  into 
damage  caused  by  overpressure  and  by  dynamic  pressure. 

The  heat  of  the  weapon  is  almost  instantly  transmitted 
through  thermal  radiation  to  every  object  within 
line-of-sight  constraints.  Objects  at  distances  of  several 
miles  will  begin  to  absorb  enough  heat  to  cause  combustable 
materials  to  burst  into  flame.  The  pressure  of  the  weapon 
is  transmitted  through  a  blast  wave  which  emanates 
spherically  from  the  point  of  detonation  and  travels  at 
speeds  closely  related  to  the  speed  of  sound.  Upon 
encountering  a  target  element  the  blast  wave  may  crush  the 
target  because  of  a  high  overpressure,  or  it  may  destroy  it 
by  translation  -  the  act  of  tumbling  it  about  on  the  ground. 

Finally,  the  nuclear  weapon  has  associated  with  it  the 
harmful,  highly  penetrating  nuclear  radiation.  The 
explosion  emits  alpha  and  beta  particles,  neutrons  and  gamma 
rays.  Due  to  their  charge,  alpha  and  beta  particles  have  no 
military  significance  in  weapon  detonation  and  are  hereafter 
ignored. 


14 


Neutrons  are  one  or  the  basic  ouildmg  blocks  of  the 
atonic  nucleus  and  have  no  electronic  charge.  It  is  the 
lack  of  charge  which  allows  deep  propagation'  and 
penetration. 

Ganna  rays  are  electromagnetic  radiation.  They  are 
identical  in  composition  to  light,  radio  waves,  and  X-rays 
except  fcr  wave  length  and  are  cot  easily  stopped  or 
absorbed. 

Beth  neutron  and  gamma  radiation  are  refered  tc  as 
ionizing  radiation.  The  biological  implication  of  this  is 
that  as  the  radiation  passes  through  living  tissue  it 
ionizes  seme  of  it,  changing  the  chemical  structure  into  a 
different,  non-f unct icnal  substance.  Essentially,  it  kills 
cells  or  inhibits  then  from  functioning  as  intended. 


B.  ENERGT  DISTRIBUTION 

The  percentages  of  total  energy  appearing  as  blast, 
thermal  radiation,  and  nuclear  radiation  depend  mostly  on 
the  altitude  at  which  the  blast  occurs.  Por  bursts  within 
the  lewer  atmosphere  the  percentages  are  about  50,  35,  and 

15  percent  respectively.  Thus,  as  with  conventional 
weapons,  blast  is  the  dominant  factor. 

C.  BLAST 

Targets  are  damaged  by  a  crushing  mechanism  caused  by 
overpressure  or  by  a  tumbling  action  caused  by  dynamic 


1  •  Cv§£^essu£s 


Overpressure  is  the  crushing  force  applied  to  a 
target.  When  an  object  begins  to  interact  with  the  blast 
wave,  the  leading  edge  of  the  object  is  subjected  to  an 
increase  in  pressure  while  the  trailing  edge  is  still  at 


ambient  pressure.  This  produces  a  quick  net  force  away  from 
the  blast.  As  the  blast  wave  envelopes  the  object  the 
entire  exterior  is  subjected  to  an  increased  pressure .  while 
the  interior  remains  at  ambient  pressure.  Thus,  a  crushing 
force  is  applied  to  the  exterior  in  an  attempt  to  damage  it. 
Lungs  and  eardrums  are  easily  damaged  by  this  mechanism  at 
overpressures  of  6  psi. 

2 .  Cynamic  Pressure 

As  the  name  implies,  this  is  the  pressure  which  is 
associated  with  high  winds.  Dynamic  pressure  can  damage 
targets  by  pushing,  tumbling  and  tearing  them  apart. 
Targets  damaged  primarily  by  dynamic  pressure  are  called 
drag  sensitive.  With  the  exception  of  heavily  armored 
vehicles  such  as  tanks,  most  military  materiel  is  drag 
sensitive.  Personnel  in  the  open  are  drag  sensitive  and  are 
damaged  easily  by  flying  debris. 


3 .  Hast  wave 


:1CS 


Opor  detonation,  a  blast  wave  of  compressed  air 
moves  spherically  away  from  the  point  of  detonation.  The 
wave  speed  is  well  in  excess  of  the  speed  of  sound.  As  the 
wave  propagates,  the  speed  diminishes  until  it  exhausts 
itself.  The  quantification  of  three  blast  wave 

characteristics  is  essential  to  simulating  blast  casualties. 
The  characteristics  are  overpressure  and  dynamic  pressure  as 
well  as  the  time  of  arrival  of  the  blast  wave. 

a.  Overpressure 


Fig.  2.1  [fief.  1:  p.  113]  is  a  graph  relating 
distance  from  ground  zero  to  height  of  burst  with 
overpressure  as  a  parameter.  The  figure  is  for  a  1  KT 
weapon  but  can  be  used  for  any  weapon  with  the  proper 
scaling  laws  applied.  Equation  2.1  is  the  cube  root  scaling 
equation  which  will  apply  here  and  in  the  next  two  sections. 


HEIGHT 


and  the  ground  distance  is 


1200 

V73 


600  ft. 


From  Fig.  2.1,  at  a  distance  of  600  ft.  with  a  height  of 
hurst  of  100  ft.  the  peak  overpressure  is  30  psi. 

Bef.  1  has  two  other  overpressure  graphs.  One 
is  fcr  overpressures  of  10,000  to  100  psi  and  the  ether  is 
for  overpressures  of  15  to  1  psi. 

These  problems  may  be  worked  in  reverse  to  solve 
for  the  optimum  height  of  burst  if  one  knows  the  required 
overpressure  to  defeat  a  target.  Phis  concept  will  be  quite 
important  in  chapters  III  and  IV.  An  example  will 
demonstrate  the  methed. 


Given:  An  8  KT  weapon  will  be  used  against  a  target  which 

is  defeated  by  10  psi  overpressure. 

Find:  Optimum  height  of  burst. 

Solution:  The  optimum  solution  is  the  one  which  produces 

the  required  effect  at  the  farthest  distance  from  the  burst. 
From  Fig.  2.1  the  distance  is  1460  ft.  and  the  height  of 
burst  is  740  ft.  From  equation  2.1  the  actual  height  of 
burst  and  actual  distance  where  the  effect  will  be  felt  can 
be  found. 

d  =  d^W1^3  =  1460*  81/3  =  2920  ft.  =  990  meters 


h  =  h1*W1/3  =  740  •81/3  =  1480  ft.  =  451  meters 

In  this  manner,  the  optimum  height  of  burst  for 
this  damage  mechanism  may  be  found  for  any  target  when  the 
required  overpressure  and  expected  yield  are  known.  This 
will  be  essential  later  in  defining  what  is  known  as  the 
governing  effect. 


t.  Dynamic  Pressure 


Fig.  2.2  (Hef.  1:  p.  117]  is  a  graph  relating 


distance  from  ground  zero  to  height  of  burst  with  dynamic 
pressure  as  a  parameter.  Like  Fig.  2.1,  it  is  for  a  1  KT 


0  ICO  200  300  400  500  600  700  800  900  1,000  1,100  1,200  1,300  t 

OISTANCE  FROM  GROUND  ZERO  (FEET) 

Figure  2.2  Scaled  Values  of  Dynamic  Pressure. 


detonation,  but  the  same  cute  root  scaling  equation  applies 
Another  example  will  demonstrate  the  method. 

Given:  a  27  KT  weapon  is  detonated  at  a  height  of  600  ft. 
Find:  The  dynamic  pressure  at  2400  ft.  from  ground  zero. 


Solution:  The  corresponding  height  of  burst  for  a  1KT  burs* 
is; 


1  '  w1/3 


and  tfce  ground  distance  is 


_  2400 


200  ft, 


=  800  ft. 


From  Fig.  2.2,  at  a  distance  of  800  ft.  with  a  height  of 
burst  of  200  ft.  the  dynamic  pressure  is  8  psi. 

In  a  manner  analogous  to  that  offered  in  the 
section  above,  an  optimum  height  of  burst  may  be  found  for 
targets  damaged  by  dynamic  pressure. 

c.  Time  of  Arrival 

Fig.  2.3  [Ref.  1:  p.  121]  is  a  graph  relating 
distance  free  ground  zero  to  height  of  burst  with  the  time 
of  arrival  of  the  blast  wave  as  a  parameter.  Like  Figs.  2.1 
and  2.2,  it  is  for  a  1  KT  detonation,  but  the  same  cube  root 
scaling  equation  applies.  Again,  an  example  will 
demonstrate  the  method. 

Given:  a  64  KT  weapon  is  detonated  at  a  height  of  2000  ft. 
Find:  The  time  of  arrival  of  the  blast  wave  at  a  distance 

of  3  miles  from  ground  zero. 

Solution:  From  previous  techniques; 


500  ft. 


d  =  3960  ft. 


From  Fig.  2.3  the  time  for  a  1  KT  burst  is  2.8  sec.  and 


1/8  1/8 

t  =  t.*Wx  =  2.8*64X/J  =  11.2  sec. 
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Figure  2.3  scaled  Values  of  1  kt  Blast  Bave  Arrival  Tiae. 

It  is  interesting  tc  note  that  as  the  yield  increases,  the 
tiae  of  arrival  increases  which  means  the  speed  of  the  blast 
wave  is  decreasing. 

D.  TBBBB1L  BADIATIOI 

Theraal  radiation  is  a  radiant  heat  transfer  mechanism 
which  propagates  with  the  speed  of  light.  Upon  burst  every 
object  with  line-of-sight  tc  the  burst  begins  to  absorb  heat 
at  a  rate  predicted  by  the  thermodynamic  laws  governing 
thermal  radiation.  Ccabustable  materials  will  ignite  if  the 
rate  and  exposure  tine  are  sufficiently  high.  Personnel 
will  suffer  first,  second,  and  even  third  degree  burns  if 
reaction  times  are  slew  and  absorption  rate  is  high.  Fig. 
2.4  [Bef.  1:  p.  291]  shows  the  relationship  between  yield, 
distance  frem  ground  zero  and  the  expected  heat  absorption. 
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Figure  2.4  Heat  Absorption  as  a  Function  of  Bange  and  Yield 
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Equation 
Nuclear 
radiation, 
destructive 
light  which 
conditions. 


2.2  is  the  linear  regression  for  Fig.  2.4. 
bursts  emit  two  distinct  pulses  of  thermal 
The  second  one  contains  about  99  percent  of  the 
power.  It  consists  of  visible  and  infrared 
is  pocrly  attenuated  by  normal  atmospheric 


Any  condition  which  affects  the  visibility  or 
transparency  of  the  air  will  modify  the  transmission  of 
thermal  radiation.  Clouds,  smoke,  fog,  rain  and  snow  will 
absorb  and  scatter  the  energy.  Also,  a  cloud  cover  above 
the  burst  can  reflect  energy  back  to  the  target  and  increase 
the  thermal  radiation  that  would  have  otherwise  traveled 
harmessly  into  the  sky.  Obviously,  any  solid  object  which 
will  not  become  a  victim  of  thermal  radiation  will  provide 
adequate  protection.  Such  objects  include,  hills,  foxholes, 
bunkers,  vehicles,  trees  and  other  people. 

As  stated  earlier,  the  exposure  time  would  be  a  factor 
in  determining  total  exposure.  However,  in  actual  tests  it 
was  disccvered  that  ccmbustable  materials  such  as  wood  frame 
houses  charred  almost  instantaneously.  Also,  due  to  the 
short  duration  of  the  thermal  pulse,  a  steady  state  thermal 
heat  transfer  condition  was  not  achieved  and  the  buildings 
were  able  tc  absorb  and  dissapate  the  initial  burst  of 
energy  and  avoid  ignition. 

The  pcint  of  the  preceeding  paragraph  is  that  reaction 
time  of  a  potential  absorber,  such  as  a  soldier  on  the 
battlefield,  will  play  an  insignificant  part  in  determining 
thermal  radiation  effects. 

Fig.  2.5  [Ref.  1:  p.  314]  shows  the  futility  of  avoiding 
burns  if  one  is  sufficiently  close  to  ground  zero.  The 
following  quote  [Ref.  1:  p.  313]  sums  it  up; 

At  the  lower  energy  yields  the  thermal  radiation  is 
emitted  in  such  a  snort  time  that  no  evasive  action  is 
possible.  At  the  higher  yields,  however,  exppsure  tc 
much  cf  the  thermal  radiation  could  be  avoided  if 
evasive  action  were  taken  within  a  fraction  of  a  second 
of  the  explosion  time.  It  must  be  remembered,  of 
course,  that  even  during  this  short  period  a  very 
considerable  amount  of  thermal  energy  will  have  been 
emitted  from  an  explosion  of  high  yield. 
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TIME  (SECONDS) 


Figure  2.5  Energy  Emmission  as  a  Function  of  Time 


B.  1QCLE1B  BA 01  AT  10  B 

1 .  Neutron  Badiation 

Upon  detonation  of  a  nuclear  weapon  a  large  quantity 
of  free  neutrons  is  released.  The  neutron  release  process 
is  ccapleted  in  less  than  a  millionth  of  a  second  [Ref.  1: 
p.  340].  This  is  where  the  term  prompt  radiation  comes 
from.  However,  due  to  collisions  in  the  early  stages  of  the 
detonation,  actual  escape  may  be  delayed  for  a  thousandth  of 
a  second  -  an  insignificant  comfort  to  potential  targets. 

Since  the  neutron  has  no  charge  it  does  not  have  a 
direct  ionizing  effect.  The  biologically  harmful  ionization 
occurs  when  the  neutron  collides  with  hydrogen  in  body 


tissue.  The  neutron  transfers  its  energy  to  the  nucleus  and 
frees  it  frca  the  accompanying  electron.  Thus  a  free 
protron  and  a  beta  particle  are  allowed  to  cause  ionizing 
tissue  damage.  Previously,  beta  had  been  disregarded  as  a 
source  of  concern.  This  was  due  to  the  impossibility  of 
getting  the  particle  into  the  body.  In  the  context  of  the 
current  discussion,  it  is  already  inside  the  body. 

Fig.  2.6  [Hef.  Is  p.  346]  is  a  graph  relating  yield 
to  the  slant  range  tc  the  burst  with  radiation  absorption  in 
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Pigure  2.6  Heutron  Radiation. 

tody  tissue  as  a  parameter.  For  a  given  yield  and  slant 
range,  the  absorbed  radiation  can  be  predicted.  Equation 
2.3  is  the  regression  equation  for  Fig.  2.6. 
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2.  Gamma  Ra dia ticn 

Upon  detonation  of  a  nuclear  weapon,  gamma  radiation 
is  emitted.  There  are  several  mechanisms  which  cause  gamma 
radiation  at  various  times  during  and  after  the  detonation. 
Ref.  1  (Sect  8.08  -  8.19)  contains  a  detailed  explanation  of 
the  mechanisms.  Rather  than  delve  into  the  sources  of  gamma 
radiation  it  is  more  important  to  predict  absorption  rates 
and  figure  out  how  to  avoid  absorption. 
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Figure  2.7  Gamma  Radiation. 


Fig.  2.7  (Ref.  Is  p.  333]  is  similar  to  Fig.  2.6. 
For  a  given  yield  and  slant  range,  the  absorbed  dose  due  to 
gamma  radiation  can  be  predicted.  Equation  2.4  is  the 
regression  equation  for  Fig.  2.7. 
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Since  gams a  radiation  is  just  light  in  a  certain 
range  c£  the  electromagnetic  spectrum,  it  behaves  in 
accordance  with  the  same  lavs,  such  as  absorption.  However, 
common  cpaque  objects  such  as  steel,  concrete,  paper,  etc. 
which  reflect  or  totally  absorb  visible  light,  only  absorb  a 
portion  of  the  gamma  radiation.  A  common  example  of  this  is 
traveling  into  a  tunnel  with  the  radio  on  in  the  car.  The 
reception  weakens  as  the  amount  of  shielding  increases. 
Table  I  lists  seme  common  materials  which  could  be  used  as 
protection. 
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TABLE  I 

Common  Gamma  Absorbers 

Material 

Tenth  Value 

Thickness  (inches) 

Steel 

3.3 

Concrete 

11 

Earth 

16 

Water 

24 

Weed 

40 

_ _ 

The  Tenth  Value  Thickness  is  that  thickness  of  the 
specified  material  which  will  cut  the  radiation  by  a  factor 
cf  10.  Thus,  3.3  inches  of  steel  or  16  inches  of  earth  will 
reduce  a  500  rad.  dose  to  50  rads.  The  addition  cf  another 
3.3  inches  cf  steel  will  further  cut  the  dose  to  5  rads. 


III.  NOCIZAB  TAB6ET  ANALYSIS 


A.  DEFINITIONS 

Before  beginning  to  make  a  weapon  selection  it  is 
necessary  to  understand  the  terns  peculiar  to  the  subject. 
Soae  of  the  more  important  ones  are  presented  here  while  the 
more  obvious  ones  are  found  in  the  glossary. 

1 .  Radius  of  Target  (R  T ) 

This  is  the  actual  radius  of  the  target  or  an 
equivalent  radius  if  the  target  is  equatable  to  a  circle. 
Ref  2  contains  a  nomograph  to  determine  the  appropriate  RT 
for  targets  equatable  to  a  circle.  More  importantly,  the 
target  elements  within  the  target  are  assumed  to  be 
uniformly  distributed. 

2.  Radius  of  Damage  (RO) 

This  is  the  distance  from  ground  zero  at  which  a 
single  target  element  has  a  50%  chance  of  receiving  the 
specified  degree  of  damage.  The  RD  is  a  function  cf  weapon 
yield,  height  of  burst,  casualty  criterion  and,  in  the  case 
of  personnel,  the  protection  level.  Fig.  3.1  shows  a  circle 
drawn  with  a  radius  from  ground  zero  such  that  as  many 
damaged  elements  (dark  dots)  are  outside  the  circle  as  there 
are  undamaged  (light  dots)  inside  the  circle.  Thus,  an 
element  at  3D  is  said  to  have  a  50%  chance  of  becoming 
damaged. 
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4 .  Circular  Error  Probable  (CEP) 


If  a  large  number  or  rounds  obeying  a  circular 
normal  distribution  were  fired  at  ground  zero  and  the  inpact 
points  plotted,  then  CEP  would  be  defined  as  that  distance 
fro®  ground  zero  that  is  exceeded  as  often  as  not.  Fig.  3.3 
shows  an  idealized  plot  of  such  firings.  &  circle  of  1  CEP 
includes  50?  of  the  rounds. 
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Figure  3.3  circular  lormal  Distribution  of  lapact. 

The  distribution  of  impacts  is  assumed  to  have  a 
bi-variate  normal  distribution  about  the  aim  point  with 
equal  variances  and  a  zero  correlation  coefficient.  Thus,  a 
relationship  between  CEP  and  variance  can  be  determined  and 


CEP  =  1.1774-a 


(3.1) 


5 .  Circular  Distrib atioa  90  (CD9 0 ) 

CD90  is  very  much  like  CEP  except  that  90%  of  the 
rounds  are  expected  to  fall  within  1  CD90  of  ground  zero. 
When  a  priority  target  is  being  considered  it  is  CD9Q  which 
is  used  as  the  dispersion  parameter.  The  relationship 
between  CD90  and  CEP  is 


CD90  =  1.83  *  CEP 


(3.2) 


6.  Probable  Errcr  Heiq  ht-of- burst  (PEH) 

Similar  to  CEE,  but  in  one  dimension,  the  PEH  is 
that  distance  above  and  below  the  desired  height  of  burst  in 
which  50%  cf  the  rounds  are  expected  to  function.  Vertical 
delivery  is  assumed  to  be  normally  distributed  abcut  the 
desired  height  of  burst.  Therefore,  the  relationship 
between  PEH  and  standard  deviation  is  easily  calculated  as; 


(3.3) 


Best  fuzing  mechanisms  are  highly  reliable  and  accurate. 
This  ensures  optimum  height  of  burst  and  maximum  casualties. 
It  also  allows  some  simplifying  assumptions  in  Chapter  IV. 

7 .  Governing  Effect 

Best  targets  are  affected  by  more  than  one  damage 
mechanism.  In  erder  to  simplify  the  analysis  the  range  at 
which  a  target  is  damaged  is  tabulated  and  the  effect 
corresponding  to  the  largest  range  is  selected  as  the 
governing  effects  Table  II  lists  hypothetical  ranges  for 
effects  for  exposed  personnel  and  personnel  ir.  open 
foxholes. 


TABLE  II 
Governing  effect 


Effect 


Overpressure 
Dynamic  pre$sui;s 
Thermal  rad3.ata.cn 
Nuclear  radiation 


Range  (meters) 


Ex|o|ed 

225 

950 

600 


Pox  hole 
575 
N  /A 
N/A 
250 


For  exposed  personnel,  it  would  appear  that  thermal 
radiation  is  the  governing  effect.  However,  due  to  the 
unpredict ibility  of  target  posture,  Thermal  radiation  is 
never  used  as  the  governing  effect  [Bef.  3:  P.  23].  Also, 
it  is  almost  universally  true  that  the  effects  of 
overpressure  are  felt  at  a  greater  range  than  the  effects  of 
dynamic  pressure.  This  leaves  only  two  criteria  to  compete 
for  governing  effect  -  blast  due  to  overpressure  and  nuclear 
radiation.  In  table  II,  the  governing  effect  for  exposed 
personnel  is  nuclear  radiation  and  for  personnel  in  foxholes 
it  is  overpressure. 

8 .  Desired  Height  of  Burst(DHOB) 

The  DHOB  is  the  maximum  of  two  possible  HOB. 

It  is  usually  desirable  to  obtain  a  low  air  burst 
without  causing  fallout.  Eg.  3.U  is  the  formula  for 
calculating  the  height  of  burst  fallout  safe  (HOBfs)  for 
yields  less  than  100  KT  and  Eg.  3.5  is  the  formula  for 
calculating  the  HOBfs  for  yields  greater  than  100  KT. 
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HOBf s  =  30  *W 


(3.4) 


HOBf s  =  55  *WJ 


H  is  the  expected  yield  plus  10*. 


(3.5) 


Obviously  this  will  produce  fallout  in  50*  of  the 
rounds.  Therefore,  a  safety  buffer  of  3.5  PEH  is  added  to 
HOBf s .  This  quantity  is  now  refered  to  as  HOB99  -  the  HOB 
at  which  there  is  a  99*  chance  of  no  fallout. 

Depending  upon  the  type  of  target  and  command 
guidance,  there  is  a  HOB  which  will  cause  optimum  damage 
(HOBopt) .  This  was  seen  earlier  in  the  discussion  on 
overpressure  and  dynamic  pressure  and  the  governing  effect 
section  above.  The  desired  HOB  (DHOB)  is  the  maximum  of 
HO 89 9  and  HOBopt.  In  most  cases,  HOBopt  is  greater  than 
HOB 9 9  and  is,  therefore,  the  listed  HOB.  It  is  only  when 
PEH  becomes  large  that  HOB99  is  greater  than  HOBopt  and  the 
DHOB  must  be  raised  above  HOBopt,  thus  reducing  expected 
coverage. 

S .  Ground  Zero  (GZ) 

The  point  on  the  ground  directly  below  the  position 
where  the  round  actually  detonates. 

10.  Desired  Ground  Zero  (DGZ) 

The  actual  aim  point.  The  DGZ  might  not  be  the 
center  of  the  target  as  will  be  seen  later.  GZ  is  only  a 
single  realization  of  the  DGZ  and  is  not  known  until  after 
the  detonation. 


B.  I AMAGE  ESTIMATION 


1 •  ££2£I£3S 

lo  cause  maximum  damage,  *.c<5  round  must  be  targeted 
at  the  center  of  the  target.  For  various  reasons,  it  may  be 
necessary  tc  deliberately  shoot  at  a  location  other  than  the 
target  center.  If  BT  and  BD  are  of  the  same  order  of 
magnitude,  then  a  less  in  coverage  will  result  as  shewn  in 
Fig.  3.4. 


A)  Less  of  coverage  B)  No  loss  of  coverage 


Figure  3.4  Loss  of  Coverage  Due  to  Offset. 

The  analyst  uses  coverage  tables  from  reference  2  in 
making  the  analysis.  Fig.  3.5  is  a  sample  of  cne  such 
table.  Each  weapon  system  and  yield  has  a  complete  set  of 
tables.  A  set  includes; 

1.  Immediate  permanent,  immediate  transient  and  latent 
lethalities  to  exposed  personnel,  personnel  in  epen  fox 
holes  and  personnel  in  tanks  (9  tables) 

2.  Moderate  damage  to  tanks,  wheeled  vehicles  and  towed 
artillery  (3  tables) 

For  each  range  and  RT  combination,  there  are  two 
coverage  figures  listed.  The  first  is  a  high  assurance 
figure  to  be  used  when  planning  ar.  attack  on  a  priority 
target.  There  is  a  90%  chance  of  obtaining  the  listed 

coverage  of  the  specified  damage  when  using  this  figure. 


The  second  figure  is  an  expected  coverage.  There  is  a  50*5 
chance  of  obtaining  the  listed  coverage  of  the  specified 
damage  when  using  this  figure. 

2.  Choosing  tfre  Correct  Bound 

It  quickly  becomes  obvious  that  there  will  be  more 
than  cne  weapon  system  which  is  adeguate.  The  question  is, 
•♦Which  cne  is  best?"  The  answer  depends  upon  current 
guidance  from  the  chain-of-ccmmand.  Several  policies  are 
possible ; 

&.  Select  maximum  yield. 

B.  Select  maximum  coverage. 

C.  Select  minimum  yield  of  the  adeguate  possibilities. 

D.  Select  minimum  coverage  of  the  adequate  possibilities. 

This  model  uses  maximum  coverage.  In  order  to 
provide  for  friendly  toop  safety,  it  will  be  seen  that 
maximum  yield  may  not  provide  maximum  coverage. 

C.  METHODS  OF  INiLYSIS 

Current  doctrine  allows  several  methods  of  nuclear 
target  analysis  ranging  from  visual  to  preclusion  oriented. 
Only  two  methods  use  detailed  mathematical  analysis.  They 
are  called  the  index  method  and  the  numerical  method. 

The  index  method  requires  that  the  desired  ground  zero 
be  the  center  of  the  target.  The  numerical  method  allows 
the  desired  ground  zero  to  be  deliberately  offset  from  the 
target  center. 

1 •  Isdex  Method 

The  index  methcd  is  the  easiest  and  quickest  method. 
To  use  this  method  the  analyst  needs  only  the  appropriate 
table  (s)  ,  the  range  from  the  applicable  weapon  delivery 
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system  and  the  target  radius.  The  appropriate  table  is 
determined  by  target  type,  desired  effects,  weapon  type  and 
weapon  yield.  Fig.  3.5  is  a  typical  table  for  personnel 
in  the  open,  latent  lethatities,  short  range  cannon,  1  KT 
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Figure  3.5  Exposed  Personnel,  Latent  Lethalities,  SRC,  1  KT. 


Entry  arguments  are  the  range  (rounded  up  to  the 
next  highest  listed  value)  and  target  radius.  If  the  target 
radius  is  less  than  the  first  listed  radius,  then  the  first 
listed  coverage  is  taken  as  the  expected  coverage.  If  the 
target  radius  is  larger  than  the  last  listed  radius  then  the 
weapon  system  is  considered  unsatisfactory.  Otherwise, 
interpolation  is  performed  on  the  radius  and  the  resulting 
coverage  is  reported  as  the  expected  initial  coverage. 

EXAKFLE.  If  a  short  range  cannon  weapon  was  employed 
against  personnel  in  the  open  with  a  target  radius  of  690 
meters  and  the  range  to  the  target  was  5200  meters  the 
expected  coverage  would  be  86.935. 
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2 .  Numerical  Method 

When  it  is  necessary  to  deliberately  shoct  -at  a 
target  with  an  aim  pcint  other  than  the  target  center,  the 
numerical  method  must  be  used.  Fig.  3.6  was  the  chart 
designed  for  usage  with  expected  co7erage,  but  has  been 

ruled  invalid  by  OSANCA.  (Fig.  3.6  may  still  be  used  for 
expected  coverage  with  no  offset.)  Therefore,  when 

employing  an  offset  aimpoiat  it  is  necessary  to  use  the  high 

assurance  chart  as  shown  in  Fig.  3.7  which  is  used  in  the 

following  discussion. 

The  numerical  method  has  the  following  6  steps. 

1.  Calculate  RD  divided  by  ST  and  locate  on  the  left  axis. 

2.  Calculate  CD90  divided  by  RT  and  locate  on  the  bottom 
axis. 

3.  The  intersection  of  the  values  predicts  the  coverage 
without  offset. 

4.  Calculate  d  divided  by  CD90  and  locate  cn  the  left 
axis. 

5.  Measure  the  distance  from  the  point  located  in  4  above 
to  the  line  labeled  " displaced  DGZ  d/CD90".  This  is 
refered  to  as  the  "measured  distance". 

6.  Apply  the  "treasured  distance"  to  the  right  of  the  point 
found  in  3  above.  This  final  point  predicts  the 
coverage  due  tc  offset.  It  must  also  be  remembered 
that  this  coverage  has  a  90%  chance  of  occurring  since 
the  solution  was  found  by  using  the  high  assurance 
table. 

The  following  example  is  depicted  in  Fig.  3.7. 
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Figure  3.6  Area  Target  Graph  (Expected  Coverage). 
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17.  SHE  HODEL 


A.  IBTBGEOCTION 

Before  the  medal  is  invoked,  a  serious  condition  must 
exist  on  the  battlefield  -  one  in  which  the  commander  with 
responsibility  and  authority  for  employment  of  nuclear 
weapons  has  decided  to  use  them.  This  decision  must  be  made 
externally  tc  the  mcdel.  Once  the  decision  is  made,  the 
model  will  make  all  of  the  decisions  necessary  to  sclve  the 
target  analysis  problem,  deliver  the  rounds  and  assess 
damages  to  target  elements. 

B.  IIPUT 

As  with  any  military  operation,  there  are  several 
essential  elements  of  information  which  must  be  available. 

1.  The  location  and  size  of  friendly  maneuver  elements 
must  be  known. 

2.  The  location  and  type  of  nuclear  capable  artillery 
batteries  along  with  their  nuclear  load  must  be  known. 

3.  The  location  size  and  type  of  enemy  targets  must  be 
estimated. 

4.  The  commander’s  guidance  for  minimum  acceptable  damage 
to  enemy  units  must  be  announced. 


C.  STEPS  OF  THE  S00E1 

Step  1.  For  each  target  select  feasible  ba ttery/yieli 
combinations . 

A.  Use  the  index  method  first. 


B.  If  the  index  method  provides  adequate  coverage, 
check  for  troop  safety.  If  offset  is  cot 
necessary,  file  as  a  solution. 

C.  If  offset  is  required,  use  numerical  method.  If 
coverage  is  adequate,  file  as  a  solution. 

D.  Of  all  feasible  solutions,  select  the  one  which 
provides  maxima  a  coverage. 

Step  2.  For  each  target,  employ  the  best  round. 

Step  3.  Assess  damage  for  each  target. 

A.  Assess  damage  to  each  tank. 

1.  Determine  the  distance  from  GZ  where  a  tank 

has  a  50 %  chance  of  damage  due  to 

overpressure . 

2.  Declare  each  tank  as  alive  or  dead  as  a 
result  of  a  Monte  Carlo  draw. 

B.  Assess  damage  for  each  troop. 

1.  Determine  the  distances  from  GZ  where  a  troop 
has  a  50*  chance  of  survival  due  to  thermal 
radiation,  overpressure  and  dynamic  pressure. 

2.  For  each  troop  in  the  target,  perforin  a  Monte 
Carlo  draw  for  each  of  the  above  damage 
mechanisms.  If  any  one  fails,  kill  the 
troop. 

3.  For  each  troop  in  the  target  which  survives 
to  this  point,  determine  nuclear  radiation 
absorption.  Place  the  trocp  in  one  cf  the  6 
categories  of; 

a.  Dead 

b.  Eermanantely  incapacitated. 

c.  Temporarily  incapacitated. 

d.  Impaired. 

e.  Latent  impaired. 


Tc  do  all  of  this  the  model  uses  various  routines  and 
functions  which  are  detailed  in  chapter  7.  The  remainder  of 
this  chapter  will  explain  how  the  model  proceeds,  what 
assumptions  are  made  and,  in  a  few  cases,  where  seme 
simplifying  approximations  are  made. 


0.  SELECT  f EASIELE  E ATTEST /YIELD  COHBIHATIOHS 

The  model  loeps  through  every  target  with  every  possible 
battery/yield  combination  in  an  attempt  to  compile  a  listing 
of  adequate  solutions.  The  proceedure  is  tc; 

1.  Use  the  index  method. 

2.  Determine  if  offset  is  necessary. 

3.  Use  numerical  method  if  offset  is  necessary 

u.  If  coverage  is  adequate,  file  as  a  feasible  solution. 

1  •  Index  Method 

Tc  start,  the  range  to  the  target  is  computed.  If 
it  is  farther  than  the  maximum  range  of  the  weapon  system 
currently  being  looped  through,  then  it  cannot  be  engaged  by 
that  battery.  The  model  will  new  loop  to  the  next  battery. 

Technically,  a  target  which  is  out  of  range  may 
still  be  attacked  successfully.  In  this  situation  an  offset 
distance  equal  to  the  difference  between  the  actual  range 
and  the  maximum  weapon  range  is  used  in  a  numerical  method 
analysis.  The  model  does  not  allow  for  such  an  occurrance. 

If  the  target  is  within  range,  the  model  loops 
through  each  yield  of  the  weapon  system.  Starting  with  the 
smallest  yield,  the  target  radius  is  compared  to  the  maximum 
table  target  radius.  If  it  is  too  large,  then  the  model 
loops  to  the  next  higher  yield.  If  it  is  not  too  large,  it 
performs  an  index  method  analysis. 


Should  the  index  method  analysis  provide  a  coverage 
which  is  egual  to  or  greater  than  the  command  guidance  for 
the  target  type,  then  a  potential  solution  has  been  found. 
If  the  coverage  is  inadequate,  then  the  model  loops  to  the 
next  higher  yield. 

It  must  be  remembered  that  the  target/battery/yield 
sequence  is  a  nested  loop  and  it  may  not  always  be  possible 
to  loop  to  the  next  higher  yield.  At  a  time  when  the 


highest  yield  for  a  weapon  is  being  looped  over,  looping  to 
the  next  higher  yield  actually  means  looping  to  the  next 
battery  and  starting  again  with  the  lowest  yield.  The  same 
applies  to  looping  to  the  next  battery  when  the  model  is  at 
the  last  battery. 

Whenever  a  feasible  solution  is  found  with  the  index 
method,  the  target  is  immediately  checked  for  necessary 
offset.  Thus  the  Icop  in  temporairly  interrupted  in  order 
to  check  for  possible  offset. 

2.  Cf fsett ing  Desired  Ground  Zero 

Recall  that  the  index  method  required  the  desired  GZ 
to  be  the  center  of  the  target.  Any  offset  necessitated  by 
troop  safety  considerations  would  require  the  use  of  the 
numerical  method  which  would  analyze  the  target  as  a  high 
priority  target  and  result  in  a  degraded  coverage.  In  the 
case  where  initial  coverage  was  adequate,  further  analysis 
must  be  dene  to  see  if  offset  is  required  and  what  effect  on 
feasability  it  has. 

Tc  perform  an  offset  is  easily  done  on  a  map  with  a 
compass,  but  is  not  as  easily  done  with  a  digital  computer. 
Fig.  4.1  shows  how  it  is  done  manually.  The  method  is  to 
draw  circles  about  each  friendly  maneuver  unit  center.  The 
radius  of  each  circle  is  the  unit  radius  plus  the  minimum 
separation  distance  (MSD)  to  insure  troop  safety.  The 
desired  ground  zero  must  lie  outside  each  circle.  The 

45 


Desired  GZ 


©  0 

Friendly  Units 


igure  4 . 1 


Manual  Solution  to  Offset  Problem 


chosen  aim  point  is  the  point  on  the  nap  which  is  closest  to 
the  target  center  and  still  exterior  to  every  circle.  - 

Simply  stated,  one  just  gets  as  close  to  the  target 
center  as  possible  while  staying  far  enough  from  each 
friendly  element  to  ensure  safety.  To  do  this  envclves 
solving  a  ncn-linear  programming  problem  with  a  non-ccnvex 
feasability  region.  The  problem  formulation  is; 

MIN  (X-X  )2  +  ( Y-Y_ ) 2 
T  T 

st.  (X-X.)2  +  (Y-Y.)2  £  (MSD  +  Unit  Radius.)2 

11  i 

i  =  1,  •  •  •  ,  No.  of  Units 

Rather  than  solve  that  problem,  it  is  easier  to 
solve  an  iterative  geometric  approximation.  In  many  cases 
the  apptcxima  ticn  is  exact.  The  approximation  is  as 
follows. 

1.  Find  the  friendly  unit  whose  perimeter  is  closest  to 
the  target  center.  If  it  is  within  MSD  of  the  target 
center,  move  the  desired  ground  zero  directly  away 
from  it  an  amount  necessary  to  satisfy  safety. 

2.  Repeat  the  process  until; 

A.  No  unit  lies  within  MSD  or 
E.  10  steps  have  been  tried. 

If  a  solution  cannot  be  found  in  10  steps  it  is 
probably  because  no  feasible  solution  exists.  Fig.  4.2  is  a 
typical  example.  Cf  the  five  units,  three  of  them  lie 
within  MSD  and  the  desired  ground  zero  must  be  moved.  Unit 
4  is  closest  and  the  desired  ground  zero  is  moved  directly 
away  from  it.  The  new  desired  ground  zero  is  acceptaole. 
It  is  also  exact  and  not  an  approximation.  This  is  the 
typical  case. 


If  unit  2  moves  forward  as  shown  in  Fig.  4.3, 


solution 


requires 


iterations  and 


approximation.  The  first  iteration  moves  the  desired  ground 
zero  tc  position  1.  However,  unit  2  is  within  MS  E  and 
desired  ground  zero  must  be  moved  again.  This  time  position 
2  is  selected.  It  is  feasible  and  close  to  the  optimum 
location. 

3 .  Numeric  al  Method 

The  numerical  method  was  discussed  in  Chapter  III, 
Section  C-2.  This  method  was  developed  for  use  by  a  nuclear 
target  analyst  in  a  field  environment  using  tools  commonly 
available  in  a  Divisional  Fire  Support  Element.  Such  tools 
are  FH  101-31-3  [Sef.  2],  a  compass,  dividers  and  a  pencil. 
Thus,  the  transposition  of  the  "measured  distance"  and 
reading  of  a  point  on  a  graph  seem  like  simple  tasks  and 
natural  things  tc  do  manually.  The  method  followed  in  the 
computer  code  is  a  slight  modification  of  the  manual  method. 

The  scales  in  Fig.  3.7  are  logrithmic.  Therefore, 
the  addition  of  the  "measured  distance"  is  really  a 
multiplication  operation  such  as; 

CD90/RT  CD90/RT  *  f  (d/CD90) 

where  f(d/CE90)  is  the  appropriate  multiplier  and  is  only  a 
function  of  d/CD90.  It  is  much  more  reliable  in  the  field 
to  allow  the  target  analyst  to  use  dividers  than  to  require 
him  tc  correctly  multiply. 

In  the  example  on  page  39,  the  initial  CD90/RT  was 
.125  and  the  displaced  CD90/BT  was  .450.  Obviously  the 
appropriate  multiplier  was  3.60. 

Fig.  3.7  has  been  included  here  as  Fig.  4.4  with 
point  number  7  added.  If  the  value  of  the  point  number  7  is 
multiplied  by  100,  it  will  be  3.60.  Thus,  the  appropriate 
multiplier  can  be  found  directly  from  rhe  figure  without  any 


further  modification.  Several  points  on  the  d/CD90  scale 
and  the  corresponding  f(d/CD90)  are  in  the  data  base.  For 
any  d/CE90  the  correct  f(d/CD90)  can  be  found  with  one 
interpolation. 

The  coverages  from  Fig.  3.7  are  input  into  the  data 
base  as  a  2-dimensional  table.  Thus  with  RD/RT  and  CD90/RT 
known,  a  double  interpolation  will  produce  the  predicted 
coverage. 

4.  Choosing  the  Best  Battery/Yield  Combination 

If,  after  all  of  this  work,  the  coverage  is  now 
inadequate,  the  solution  is  abandoned  and  the  program  moves 
cn.  However,  if  it  is  adequate,  a  solution  containing  the 
target,  battery,  yield,  desired  SZ  and  predicted  coverage 
can  be  generated  and  filed  away  for  later  use.  It  would  be 
useful  to  file  the  feasible  solution  in  such  a  way  that 
recovery  cf  the  top  one  on  the  list  is  the  best  solution. 

Table  III  shews  hypothetical  results  of  an  analysis 
requiring  offset  for  the  rounds  of  the  MRC  battery.  Dus  to 
safety  imposed  offset,  the  largest  yield,  while  still 
meeting  command  guidance,  is  clearly  not  the  best  choice. 
In  fact,  neither  of  the  rounds  from  the  MRC  battery  is  the 
best  choice.  The  SRC  battery  with  a  1  KT  yield  will  provide 
the  highest  coverage.  Therefore,  this  program  will  file 
solutions  in  descending  order  cf  coverage.  Should  a  tie 
occur,  the  SRC  takes  precedence  over  the  MRC  for 
conservation  of  force  reasons. 

£.  EMPLOYING  THE  ROQHD 

When  all  of  the  targets  have  been  analyzed  for  each 
battery,  the  target  analysis  portion  of  the  model  is 
finished.  It  is  time  to  select  the  best  battery/yield 
combination,  employ  the  weapon  and  assess  damage.  Actual 


TABLE  III 

Degraded  Coverage  Due  to 

Offset 

Battery 

Type 

Initial 

coverage 

Offset 
cov  erage 

Short  Eange 
Cannon 

.2  Kt 
1  KT 

41% 

78% 

4  1% 

76% 

Hedium  Range 
Cannon 

2  KT 
8  KT 

82% 

97% 

•p 

cu 

MM 

1 

1 

burst  parameters  are  normally  distributed.  The  actual 
ground  zero  is  picked  from  a  bivariate  normal  distribution. 

G2  =  N  (desired  ground  zero,  CEP/1.1774) 

Hieght  cf  burst  and  yield  are  selected  from  univariate 
normal  distributions. 

HOB  *  N  (desired  HOB,  PEH/.67) 

Yield  =  N  (nominal  yield,  nominal  yield/10) 


F.  ASSESSING  DAHAGE 

For  damage  assessment  it  is  assumed  that  HOBopt  was 
selected  as  though  overpressure  was  the  governing  effect. 
This  is  dene  to  avoid  reading  Figs.  2.1  and  2.2  into  the 
data  base.  This  is  dene  for  two  reasons. 

1.  The  data  in  the  figures  is  very  coarse  and  should  not 
be  accepted  as  highly  accurate  [P  f.  1:  preface]. 

2.  A  linear  function  results  from  the  assumption. 


If  the  assumption  that  the  weapon  will  function  at  HCBcpt 
for  ewer  pressure  is  made,  then  one  variable,  the  HOB,  can  be 
omitted.  If  the  employed  weapon  functions  at  HOBopt  then 
the  scaled  HOB  is  sixply  the  HOBopt  for  a  1  KT  weapon. 

Using  this  HOB  and  the  ranges  found  in  Pig.  2.  1  for 
desired  overpressures  and  in  Fig.  2.2  for  desired  dynamic 
pressures,  it  is  possible  to  plot  yield  against  range  with 
pressure  as  a  parameter.  Figs.  4.5  and  4.6  show  these 
results. 
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Figure  4.5  Overpressure  from  Optimum  HOB. 


If  the  weapon  yield  and  the  required  overpressure  or 
dynamic  pressure  for  a  particular  target  defeat  is  known,  it 
is  a  simple  matter  for  rhe  computer  to  calculate  the  range 
at  which  defeat  will  occur.  An  example  will  demonstrate 
this. 
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Figure  4.6 


Yield  (KT) 

Dynamic  Pressure  from  Optiaua  HOB. 


Given:  A  detonation  results  in  a  10  KT  yield.  An  AFC  is 

damaged  fcy  15  psi  overpressure. 

Find:  The  distance  from  ground  zero  at  which  the  AFC  is 

damaged. 

Solution:  From  Fig.  4.5  locate  10  KT  on  the  lower  scale. 

Go  up  to  the  15  psi  line.  Bead  790  meters  from  the  left 
scale. 

In  the  model.  Eg.  4.1  is  used  to  determine  the  range. 
The  coefficient  A  is  a  function  of  overpressure  or  dynamic 
pressure  and  is  determined  from  Figs.  4.5  and  4.6  when  the 
yield  is  1  KT. 


A* Yield' 


(4.1) 


The  same  simplification  can  be  made  concerning  arrival 
time.  Equation  4.2  is  the  result. 


T  =  { ( .  0028  *R)  -  .45} *Yield1/3 


(4.2) 


It  must  be  pointed  out  that  this  method  will  predict  a 
larger  8D  for  overpressure  than  is  correct  when  overpressure 
is  net  the  governing  effect.  However,  overpressure  is 
almost  always  the  governing  effect.  (Radiation  is  usually 
the  governing  effect  only  against  personnel  in  the  open.) 
But,  when  radiation  is  the  governing  effect,  enhancing  the 
effects  of  overpressure  to  their  optimum  will  still  not 
outweigh  the  effects  of  radiation.  k  soldier  killed  by 
radiation  will  still  be  killed  by  radiation.  The  only  error 
is  that  a  soldier  who  survives  radiation  and  should  have 
survived  overpressure  will  now  have  a  decreased  chance  of 
survival. 

The  alternative  to  the  simplified  method  is  to  revert  to 
Eg.  2.1  with  a  table  look-up  and  2  way  interpolation. 


G.  71BI1BIIITT  OP  EFFECTS 

An  object  positioned  at  a  certain  range  from  a  given 
yield  will  be  subjected  to  overpressure,  dynamic  pressure, 
thermal  radiation  and  nuclear  radiation.  In  many  cases  one 
damage  mechanism  will  dominate  the  others.  It  would  be  a 
simple  matter  to  use  only  the  dominant  mechanism  in 
determining  damage.  In  an  expected  value  model  this 
technique  is  entirely  justified  but  such  a  proceedure  is  not 
followed  in  this  model. 

In  a  probabilistic  model,  such  as  this  one  it  is 
entirely  possible  for  a  target  element  to  "beat  the  odds"  on 
overpressure  and  be  killed  by  secondary  missiles  brought 
about  by  dynamic  pressure.  If  there  was  a  75*  chance  of 


damage  by  overpressure  and  only  a  10*  chance  of  damage  by 
dynamic  pressure  in  a  50*  threshold  model,  the  element  would 
be  killed  by  overpressure  every  time.  In  this  model,  each 
mechanism  gets  a  Monte  Carlo  chance  to  cause  damage.  In  the 
above  example  the  subject  has  a  77.5*  chance  of  damage. 
More  importantly,  it  has  a  non-zero  chance  of  survival. 

If  the  HD  from  Fig.  4.5  is  found  to  be  700  meters  (10 
psi  and  10  KT)  then  what  are  the  probabilities  for  survival 
at  a  closer  range  and  for  damage  at  a  fartner  range?  The 
expected  value  model  would  set  these  probabilities  at  zero. 
This  model  will  allow  survival  at  a  closer  range  and  damage 
at  a  farther  range.  To  do  this  requires  a  probability 
distribution. 
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Figure  4.7  Variability  of  Effects. 


Fig.  4.7  [Ref.  4:  p.  D-3]  is  used  in  the  FM  101-31 
series  and  is  followed  in  this  model.  The  20*  variance 
curve  is  the  one  which  is  used.  Plotted  on  normal  clotting 
paper,  it  locks  like  Fig.  4.8. 


According  to  Fig.  4.7 ,  a  target  element  at  ground  zero 
has  some  finite  chance  of  survival  while  an  element  at  an 
infinite  range  may  still  perish.  This  is  a  logical  result 
of  the  normal  probability  assumption.  To  approximate 
reality,  any  element  with  less  than  0.5%  chance  of  survival 
is  killed  while  anything  with  more  than  99.5%  chance  of 
survival  always  survives.  Thus,  only  when  R/RD  is  between 
.5  and  1.5  is  the  Monte  Carlo  chance  taken.  An  example  will 
demonstrate  how  the  method  works. 

GIVES:  A  6  KT  burst.  Tanks  are  damaged  by  20  psi 

overpressure . 

FIND:  Damage  to  tanks  at  200  and  900  meters  from  GZ. 

SOLUTION:  From  Fig.  4.5  the  range  at  which  a  tank  is 

damaged  is  500  meters.  Correctly  stated,  this  is  RD,  the 
distance  at  which  the  tank  has  a  50%  chance  of  survival  The 
first  tank  has  a  distance/RD  of  .40.  From  Fig.  4.3,  the 

probability  of  survival  is  less  than  0.5%  and  the  tank  is 

killed. 

The  ssccnd  tank  has  a  distance/RD  of  .555.  A  uniform 

(0,1)  draw  is  made  for  the  tank  and  the  number  of  .60  is 

drawn.  Since  .555  is  less  than  .60,  the  tank  survives.  Had 
.554  cr  less  bean  drawn,  the  tank  would  have  been  killed. 

H.  TARGET  RESPONSES 

Target  elements  are  of  two  types  -  tanks  and  personnel. 
However,  any  target  element  which  is  damaged  by  blast, 
thermal  radiation  or  nuclear  radiation  may  be  easily  input 
into  the  model.  Such  elements  would  include  APC's,  trucks, 
buildings,  aircraft,  forests,  bridges,  etc. 

In  this  model  the  tanks  are  either  fully  functional  or 
killed.  Personnel  may  have  an  additional  state  -  dying. 
This  is  caused  by  the  absorption  of  nuclear  radiation. 
While  the  various  mechanisms  may  kill  a  soldier  upon  burst. 
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only  the  nuclear  radiation  will  causa  a  latent,  lingering 
death.  Fig.  4.9  [Hef.  4:  p.  C-9]  is  used  to  classify 
radiation  casualties  into  the  states  of: 

1.  Dead. 

2.  Incapacitated  and  dying. 

3.  Iapaired  and  dying. 

4.  Safe  and  fully  functional. 
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Figure  4.9  Eersonnel  fiesponses  to  fiadiation. 

The  model  simplifies  Fig.  4.9  into  Fig.  4.10  wi^h  the 

following : 

1.  Anyone  receiving  18,000  rads  is  declared  dead 
immediately. 

2.  Anyone  receiving  8,000  to  18,000  rads  is  immediately 
and  permanantly  incapacitated.  He  cannot  shoot,  move  or 
communicate  but  is  still  alive  and  should  be  considered 
a  potential  target. 

3.  Anyone  receiving  3,000  to  3,000  rads  is  immediately 
incapacitated  but  in  time  is  upgraded  to  functionally 
impaired.  The  time  is  uniformly  distributed  from  30  to 
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45  minutes.  Eventually  the  person  dies.  Death  occurs 
in  a  uniformly  distributed  amount  of  time  from  2  to  5 
days. 

4.  Anyone  receiving  2,000  to  3,000  rads  is  immediately 
placed  in  the  functionally  impaired  state  where  he 
remains  until  death  occurs  in  2  to  5  days. 

5.  Anyone  receiving  a  dose  greater  than  his  own  personal 
lethal  limit  bet  less  than  2,000  rads  is  unaffected 
until  some  latent  period  has  passed  at  which  time  he 
becomes  functionally  impaired.  This  time  is  uniformly 
distributed  from  30  to  60  minutes.  Eventually,  the 
person  dies  in  2  to  5  days.  The  lethal  limit  is 
normally  distributed  about  650  rads  with  a  standard 
deviation  of  50  rads. 

6.  Anyone  receiving  less  than  his  own  personal  lethal 
limit  is  unaffected. 

I.  OUTPUT 

When  the  model  is  completely  finished  in  will  print  all 
of  the  pertinent  information  concerning  each  target  and 
employed  round.  Final  output  will  be; 

1.  Burst  parameters  of  GZ,  HOB  and  yield  for  destruction 
of  each  target. 

2.  Location  and  disposition  of  each  element  in  each 
target.  For  personnel  not  immediately  killed,  the 
incapacitation  and/or  impairment  times  will  be  listed  as 
appropriate. 
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Figure  4.10  Simplified  Personnel  Responses  to  Radiation. 


J.  SUMH1BY 

This  model  will,  given  several  nuclear  capable  artillery 
batteries,  friendly  locations  and  enemy  locations,  select 
optimum  choices  of  batteries  and  yields  to  fire  at  each 
enemy  target.  It  will,  when  necessary,  allow  fcr  troop 
safety  by  offsetting  the  aim  point  of  the  weapon.  It  will 
simulate  delivery  of  each  round  and  stochastically  assess 
damage  to  each  individual  element  within  the  target.  Tanks 


V.  MODEL  BOPTIHES 


This  chapter  will  explain  the  important  sections  of  the 
actual  model  in  detail.  Actual  lines  of  the  SIHSCRIPT  model 
(printed  in  capital  letters)  will  be  presented  followed  by 
an  explanation  of  what  they  do  or  mean.  The  complete  modal 
is  included  as  Appendix  B.  Each  routine  and  function  is 
discussed  as  a  sepatate  section  in  this  chapter. 

In  reading  about  each  routine  and  function  it  may  be 
helpful  tc  refer  to  Fig.  5.1  to  keep  the  flow  of  the  model 
in  order. 

A.  PREAMBLE 

The  preamble  extablishes  global  variables,  identifies 
permanent  and  temperary  entities  and  their  attributes, 
groups  entities  into  sets  and  identifies  user  defined 
functions. 

PERMANENT  ENTITIES 

EVERY  EATTEBY  HAS  A  SIZE,  AN  X3,  A  Y3  AND  SOME  NOC. ROUNDS. 

EVERY  COMPANY  HAS  A  XC ,  A  YC,  A  ZC  AND  A  CO. RADIUS 

EVERY  TARGET  HAS  A  TYPE,  A  XT,  A  YT,  A  ZT,  A  P.T 

AND  MAY  OWN  A  TANK.  SET,  A  TROOP. SET,  AND  A  LISTING 

The  EATTERY  is  an  entity  representing  a  nuclear  capable 
artillery  battery.  The  SIZE  is  either  short  range  cannon 
(155mm,  SIZE  =  1)  or  medium  range  cannon  (8  in.,  SIZE  =  2). 
XB  and  YE  are  grid  coordinates.  NOC. ROUNDS  is  the  quantity 
cf  nuclear  rounds  the  BATTERY  has  available  to  fire. 
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Figure  5.1  Routine/Function  Routing  chart. 

The  COMPANY  is  an  entity  representing  a  friendly  unit. 
It  may  be  an  individual  soldier  or  vehicle  or  a  collection 
cf  such.  Its  reason  for  existing  is  to  cause  offset 
targeting  due  to  friendly  safety  constraints.  xc,  YCr  ZC 
are  the  grid  coordinates  and  altitude.  CO. RADIUS  is  the 


radius  of  the  unit. 


The  TARGET  is  an  entity  which  is  to  be  attacked.  XT,  YT 
and  ZT  are  grid  coordinates  and  altitude.  RT  is  the  target 
radius  cr  equivalent  as  explained  in  Chapter  III.  The 
TARGET  owns  three  sets.  They  are  sets  of  TANKS,  TROOPS  and 
a  LISTING  of  feasible  SOLUTIONS. 


TEMPORARY  ENTITIES 

EVERY  TANK  HAS  A  X  .TANK,  A  Y . TANK,  A  TANK. DEAD, 

AND  HAY  BELONG  TO  A  TANK.  SET 
EVERY  TROOP  HAS  A  X. TROOP,  A  Y. TROOP,  A  TRA NS . FACTOR , 
A  TROOP. DEAD,  A  LETHAL. DOSE,  AN  ACCOM. DOSE, 

A  T.LETH,  A  1. IMPAIR. NOC 
AND  MAI  BELONG  TO  A  TROOP. SET 
EVERY  SOLUTION  HAS  A  UNIT .TO. FIRE,  AN  XDGZ,  A  YDGZ, 

AN  I  YIELD,  A  JROW,  A  KCCL  AND  A  PCT 
AND  BELONGS  TO  A  LISTING 


The  TANK  has  coordinates  X. TANK  and  Y . TANK.  The 

altitude  is  assumed  equal  to  that  of  the  target  to  which  it 
will  be  assigned.  TANK. DEAD  is  an  integer  number,  either  0 
cr  3.  AO  means  the  tank  is  fully  functional.  A  3  means  it 
has  suffered  a  catastrophic  kill. 

The  TROOP  has  coordinates  X. TROOP  and  Y. TROOP.  The 

altitude  is  assumed  equal  to  that  of  the  target  to  which  it 
will  be  assigned.  The  TRANS. FACTOR  is  the  radiation 
transmission  factor  and  is  currently  either  a  1  for  exposed, 
.3  for  protection  inside  an  APC  and  .00  1  for  a  fcx  hole. 
TROOP. DEAE  is  similar  to  TANK. DEAD  except  for  two  additional 
values.  A  1  is  functionally  impaired  and  a  2  is 

incapacitated.  LETHAL.  DOSE  is  that  radiation  dose  above 
which  the  TROOP  will  eventually  die.  ACCUM.DOSE  is  the 
radiaticn  dose  the  TECOP  has  currently  accumulated.  T.LETH 
is  the  time  when  TROOP. DEAD  is  to  be  set  equal  tc  3. 


T. IMPAIR. NUC  is  the  time  when  TROOP. DEAD  is  to  be  set  equal 
to  1 . 

A  SCIUTION  is  an  answer  to  the  nuclear  target  analysis 
problem.  Its  attributes  are  the  firing  unit,  desired  ground 
zero,  expected  coverage  and  indicies  to  locate  the  table, 
row  and  column  where  the  solution  was  found.  It  is  filed  in 
a  LISTING  belonging  to  a  TARGET.  The  LISTING  is  ranked 
first  by  high  coverage  and  then  by  low  yield. 

GLOBAL  VARIABLES 

There  are  two  sets  of  equivalent  arrays.  one  pertains 
to  the  SRC  and  the  ether  to  the  MRC.  The  SRC  has  9  rows 
while  the  MRC  has  19  rows.  In  the  following,  only  the  SRC 
is  explained  but  the  description  extends  equally  to  the  MRC. 

SRC. E .COVER AGE  is  the  4  dimensional  matrix  which 
contains  the  coverage  probabilities.  The  first  subscript  is 
the  yield.  The  second  is  the  type  of  target.  The  third  and 
fourth  are  the  row  and  column  of  the  table  described  by 
yield  and  type.  The  current  model  allows  seven  types  of 
targets.  They  are; 

1.  Exposed  personnel.  Immediate  permanent  casualties. 

2.  Expcsed  personnel.  Immediate  transient  casualties. 

3.  Personnel  in  fox  holes.  Immediate  permanent  casualties. 

4.  Personnel  in  fox  holes.  Immediate  transient  casualties. 

5.  Moderate  damage  to  Tanks. 

6.  Moderate  damage  to  Wheeled  Vehicles. 

7.  Moderate  damage  to  Artillery. 

SRC. RANGES  is  the  array  of  ranges  that  the  system  is 
capable  cf.  (See  Fig.  3.5) 

SRC. TABLE. RT  is  the  3-d imsnsiona 1  matrix  which  contains 
the  RT  fer  the  various  tables.  The  first  subscript  is  the 
yield.  The  second  is  the  type  of  target.  The  third  is  the 
column  of  the  table.  (See  Fig  3.5) 


SBC. YIELDS  is  the  array  of  possible  yields.  Each  weapon 
has  two  yields  available.  The  SBC  has  a  . 2  KT  and  a  1KT 
while  the  MfiC  has  a  2  KT  and  an  8  KT. 

USD. SBC  is  a  2-dimensional  matrix  which  contains  values 
of  minimum  safe  distances.  The  first  subscript  is  the  row 
of  the  table  and  the  second  is  the  yield.  MSD  is  cnly  a 
function  cf  weapon  size,  range  and  yield  but  not  target  type 
or  degree  of  damage. 

CD90.SRC,  CEP. SBC  and  PEH.SRC  are  arrays  with  only  the 
range  as  a  variable.  They  are  a  function  cf  the  delivery 
system  cnly.  The  subscript  is  the  row  of  any  table  for  the 
applicable  weapon  system. 

HCB.SBC  is  a  3-dimensional  matrix  which  has  the  desired 
height  of  burst.  Tte  first  subscript  is  the  yield.  The 
second  subscript  is  the  row  of  the  table.  The  third 
subscript  is  the  type  of  target. 

SBC. MIN. RD  and  SBC. EXP. HD  are  the  minimum  radius  of 
damage  and  the  expected  radius  of  damage.  Their 

dimensicnality  and  subscript  meaning  are  the  same  as 
HOB. SBC. 

CD90.RT  ,  RD.RT,  D. CD9 0  and  DISPLACED. COVERAGE  are  the 
parameters  used  with  Fig.  3.7  to  solve  an  offset  coverage 
problem.  CD90.RT  is  CD90/RT.  RD.RT  is  RD/RT.  D.CD90  is 
d/CD9Q.  DISPLACED  COVERAGE  is  a  matrix  of  coverage  values 
used  for  table  look-up  once  the  entry  arguments  have  been 
determined. 

B.  MAIN 

The  MAIN  serves  as  a  vehicle  to  define  and  read  local 
variables  and  call  subroutines  and  functions  which  dc  the 
bulk  cf  the  work. 
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POLICY  is  an  array  which  represents  the  commander* s 
guidance  cn  each  of  the  types  of  targets.  Its  length  is 
NTABLES  -  one  for  each  of  the  seven  types  of  targets.  If 
the  predicted  coverage  on  a  target  does  not  equal  or  exceed 
the  POLICY  for  that  type  of  target,  then  the  solution  is 
abandcned. 

RANGE  is  the  range  from  the  battery  to  the  target  for 
which  a  solution  is  being  sought. 

The  remainder  cf  MAIN  is  involved  with  calling 
subroutines  to  solve  the  problem.  MAIN  first  calls 
READ. RESERVE  which  reserves  and  reads  all  of  the  global 
variables.  MAIN  then  calls  INITIALIZE.  It  then  loops 

through  every  TARGET  and  every  BATTERY  in  an  attempt  to  find 
adequate  SOLUTIONS.  To  do  this  it  calls  ANALYSIS  which,  in 
turn,  calls  other  subprograms.  After  returning  from 
ANALYSIS,  the  MAIN  is  ready  to  employ  the  weapons.  It 
removes  the  first  SOLUTION  from  the  LISTING  of  each  TARGET 
and  calls  DETONATION  which  assesses  damage.  Upon  returning 
from  DETONATION,  the  MAIN  lists  the  attributes  of  the 
permanent  entities  and  terminates. 

C.  SUBROUTINE  READ.  RESERVE 

This  subroutine  reserves  and  reads  the  global  variables 
as  defined  in  the  PREAMBLE 

D.  SUBROUTINE  INITIALIZE 

This  is  a  surrogate  for  a  driver  program.  As  stated  in 
the  abstract,  this  program  is  intended  as  a  subroutine  for  a 
larger  simulation  -  cne  which  already  has  entities  such  as 
TARGET,  TANK  and  TROOP.  Since  these  entities  do  not  exist, 
INITIALIZE  creates  them.  A  detailed  description  of  the 
subroutine  will  not  be  offered  here  since  it  is  unimportant 
and  easy  to  read.  It  will  suffice  to  say  that  if  uniformly 
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distributes  TANK  and  TROOP  entities  within  the  RT  of  the 


TARGET  while  assigning  attributes  to  the  temporary  entities. 

E.  SUBROUTINE  ANALYSIS 

This  subroutine  is  the  origin  of  SOLUTION.  It  is  here 
that  a  decision  is  made  to  create  a  SOLUTION  and  then  to 
file  it  in  a  LISTING  or  discard  it.  The  first  acticn  is  to 
determine  the  size  cf  the  EATTERY  which  has  been  passed  to 
it.  It  then  determines  if  the  TARGET  is  within  range.  If 
not,  it  returns  to  BAIN  without  a  SOLUTION.  If  the  range  is 
less  than  the  maximum  weapon  range,  then  it  finds  the  proper 
row  (J)  to  use  as  an  entry  argument.  It  then  loops  ever  all 
yields  tc  determine  coverage  for  each  yield.  For  each  yield 
it  determines  the  proper  column  (K)  to  use  as  the  ether 
entry  argument.  If  it  cannot  find  any  (RT  too  big)  if  loops 
to  the  next  yield  cr  returns  to  BAIN  if  it  was  at  the 
largest  yield. 

If  a  row  and  column  were  found,  then  a  SOLUTION  is 
created  and  the  index  method  is  used  to  compute  the  initial 
coverage.  If  the  ET  is  less  than  or  equal  to  the  first 
listed  in  the  table,  then  the  coverage  in  the  first  column 
is  reported.  If  not,  the  function  EXP. COVERAGE  is  called. 
This  returns  the  expected  initial  coverage  using  the  index 
method. 

A  decision  is  made  concerning  the  coverage  and  the 
command  guidance  (POLICY)  for  that  type  of  target.  If  the 
coverage  is  adequate,  then  OFFSET. COVERAGE  is  called.  This 
will  return  with  a  new  coverage  which  is  less  than  cr  equal 
to  the  cne  returned  by  EXP  .COVERAGE.  If  this  new  coverage 
is  still  adequate,  the  attributes  of  SOLUTION  are  determined 
and  it  is  filed  in  the  LISTING  of  the  TARGET  for  which  if 
was  called.  Control  is  then  returned  to  MAIN. 
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F.  SOBBOOTIME  EXP.  COVES  AGE 


This  function  was  originally  established  to  perforin  the 
interpolation  with  the  index  method.  However,  as  it  became 
necessary  tc  perform  interpolation  in  other  subprograms,  a 
dedicated  interpolation  function  (INTEBP)  was  written. 
Currently  this  function  simply  determines  the  size  of  the 
BATTERY  and  properly  formats  the  arguments  for  a  call  to 
INTERE.  It  then  returns  with  the  expected  coverage  to 
ANALYSIS. 

G.  SOBHOOTINE  OFFSET .COTES AGE 

This  subroutine  is  called  when  the  initial  analysis  is 
at  least  adequate.  It  determines  if  any  offset  is  necessary 
and,  if  so,  performs  an  analysis  using  the  numerical  method. 

The  first  order  of  business  is  to  call  MAX. DISTANCE 
which  actually  locates  the  new  desired  ground  zero.  If 
there  is  no  displacement,  then  there  is  no  change  in 
coverage  and  control  is  returned  to  ANALYSIS.  max. distance 
will  also  return  a  flag  which  may  mean  no  feasible  solution 
exists  and  coverage  is  set  equal  to  zero.  This  will  cause 
the  solution  not  to  be  filed  in  a  LISTING. 

Next,  the  battery  size  is  determined  and  the  entry 
arguments  of  Fig.  3.7  are  calculated.  Depending  on  the 
value  of  d/CB90  (Y2)  one  of  three  things  can  happen; 

1.  Consider  the  target  as  a  point  target. 

2.  Analyze  as  instructed. 

3.  Degradation  is  insignificant.  Return. 

Currently,  a  point  target  analytical  technique  is  not 
included,  tut  the  technique  is  slightly  easier  than  this 
one.  The  only  time  an  area  target  would  be  analyzed  as  a 
point  target  is  if  the  RD  is  at  least  10  times  greater  than 
RT.  This  only  happens  with  large  yields.  Case  2  is  the 


7  1 


required  (assuming 


only  one  for  which  further  analysis  is 
case  1  dees  not  happen) . 

Case  2.  d/CD90  is  located  on  the  ordinate  and  (with  a 
call  tc  INTERP)  the  measured  distance  is  found.  As 
explained  in  Chapter  III,  the  measured  distance  (DELX)  is 
multiplied  by  the  original  CD90/RT  and  used  to  find  the 
predicted  coverage  due  to  offset.  A  double  interpolation  is 
necessary  to  find  the  coverage.  Control  is  then  returned  to 
ANALYSIS.  If  the  coverage  is  still  adequate,  the  SOLUTION 
is  filed  in  a  LISTING. 

H.  FUNCTION  INTERP 

This  is  an  elementary  interpolation  function,  is  self 
explanatory  and  is  net  further  explained  here. 

I.  SUBROUTINE  H  AX.  DISTANCE 

This  subroutine  performs  the  approximation  to  the 
non-linear  programming  problem  as  stated  in  Chapter  III. 
The  logic  fellows  the  graphical  solutions  as  shown  in  Figs. 
4.1  and  4.2 

J.  SUBROUTINE  DETONATION 

The  purpose  of  this  subroutine  is  to  assess  damage  now 
that  the  analysis  is  complete.  The  best  SOLUTION  for  a 
TARGET  has  teen  removed  from  its  LISTING  and  passed  zo  rhis 
subroutine. 

Again,  a  decision  is  made  regarding  the  BATTERY  size. 
The  standard  deviation  for  ground  zero  location  is  found  Eq 
3.1.  The  PEH  and  desired  HOB  are  selected  from  the  tables. 

To  determine  ground  zero  the  radial  miss  distance,  R,  is 
computed  as  N(0 , CEP/1. 1774)  while  an  angle.  A,  is  picked 
from  0(0,2  f  ).  Ground  zero  is  then  placed  at  the  polar 
coordinates  of  (R,A) 


The  actual  yield  is  selected  from  a  normal  distribution 
with  nominal  yield  as  the  mean.  (Actual  variances  are 
classified)  The  standard  deviation  is  assumed  to  be  one 
tenth  of  the  mean. 

HOB  is  selected  from  a  normal  distribution  vith  DHCB  as 
the  mean  and  PEH/.67  as  the  standard  deviation. 

The  actual  GZ,  yield  and  HOB  have  nov  been  fixed  and  the 
burst  parameters  are  output  by  the  statement; 

LIST  1,  Y ,  YIELD,  HOB 

EACH  TANK  in  TANK .  SET  (TARGET)  and  TROOP  in 

TROOP. SET  (TARGET)  is  examined  to  check  for  damage.  As 
stated,  the  TANK  can  only  be  fully  operatinal  (TANK. DEAD  = 
0)  or  killed  (TANK.  DEAD  *3).  The  only  mechanism  which  has 
any  reasonable  chance  to  destroy  a  tank  is  overpressure. 
The  overpressure  at  which  a  TANK  has  a  50%  chance  of 
survival  is  25  psi.  Thus; 

LET  RD. OVER  .PRESSURE  *  OVER  .  PRESSURE  (YIELD,  25) 

will  return  the  distance  from  GZ  where  25  psi  will  be 
encountered  and  will  be  used  to  determine  if  the  TANK  is 
destroyed. 

IF  (R/RD. OVER.  PRESSURE  LE  .5)  OR 

(.01*  (N20TRCN(R  ,  YIELD)  PLUS  GAMMA  (R,  YIELD)  )  GT  18000) 

is  a  compound  IF  statement  which  destroys  the  TANK 
(TANK. DEAD  =  3)  if  either  of  the  following  conditions  are 
met . 

1.  If  the  probability  of  destruction  by  overpressure  is 
greater  than  99.5%  or 

2.  If  the  crew  inside  would  receive  more  than  18,000 
rads  of  radiation. 
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The  TBOOP  is  put  through  all  four  damage  mechanisms 
before  being  declared  a  survivor.  The  thermal  level  at 
which  a  TROOP  has  a  50%  chance  of  survival  is  20  c al/cm  . 
Thus ; 

LET  RD. RE AT  *  REAT(yield#20) 

will  return  the  distance  from  GZ  where  20  cal/cm  will  be 
encountered  and  will  be  used  to  determine  if  the  TROOP  is 
killed.  Tc  become  a  thermal  casualty  a  compound  IF 
statement  must  be  true: 

The  TROOP  must  be  exposed  (TRANS. FACTOR  (TROOP)  *  1)  and 
The  TRCCP  must  receive  enough  thermal  energy  to  become  a 
therxal  casualty  (PROB  ABILITY(R/RD  .HEAT)  LT  CHANCE)  and 
His  reaction  time  must  be  greater  than  .  1  second. 

If  the  TROOP  survives  the  first  test,  he  may  beccme  a 
casualty  due  to  overpressure.  The  overpressure  at  which  a 
TROOP  has  a  50%  chance  of  survival  is  10  psi.  Thus; 

LET  RC. OVER. PRES  SORE  *  OVER .PRESSURE (YIELD , 10) 

will  return  the  distance  from  GZ  where  10  psi  will  be 
encountered 

and  will  be  used  to  determine  if  the  TROOP  is  destroyed. 

If  R/BD  is  less  than  .5,  the  TROOP  is  killed  (TROOP. DEAD  * 
3).  Otherwise  he  must  take  a  Monte  Carlo  chance; 

IF  PROBABILITY  (R/RD  .  OVE  R.  P  RESSURE)  LT  CHANCE 

If  the  above  statement  is  true  then  the  TROOP  is  killed. 

If  the  TROOP  is  still  alive,  he  must  now  survive  the 
effects  cf  dynamic  pressure.  This  section  is  identical  tc 
overpressure  with  the  exception  of  a  call  to 
DYNAMIC. PRESSURE  instead  Of  OVER. PRESSURE  to  get 
RD.  DYNAMIC.  PRESSURE 


After  all  of  this,  if  the  TROOP  is  still  alive,  nuclear 
radiation  gets  its  chance.  The  combined  effects  of  neutron 
and  gamma  gamma  radiation  are  added  to  the  current  radiation 
level. 

ADD  TRANS. FACTOR(TRCOP)  *((  NEUTRON  (R  ,  YIELD)  +GAMMA  (R,  YIELD))  ) 
TO  ACCOM.  DOSE(TROOP) 

If  the  IROCP  has  accumulated  more  than  18,000  rads  he  is 
killed  (TROOP.  DEAD  *  3). 

If  the  ACCOM. DOSE  is  between  8,000  and  18,000  rads  he  is 
placed  in  the  incapacitated  state  (TROOP.  DEAD(TROOP)  =  2) 

until  he  is  killed  in  2  to  5  days. 

LET  T.  LETH  (TROOP)  =  RANDI .F  (  2*24*60, 5 *24*60 , 1) 

If  the  ACCOM. DOSE  is  between  3,000  and  8,000  rads  he  is 
placed  in  the  incapacitated  state  for  30  to  45  minutes. 

LET  T. IMPAIR. NOC  (TROOP)  =  R  ANDI.  F  (30 , 45,  1) 

At  the  end  of  this  time  he  is  placed  in  the  functionally 

impaired  state  (TROOP. DE AD ( TROOP)  *  1)  .  He  is  killed  in  2 
to  5  days. 

If  the  ACCOM. DOSE  is  between  2,000  and  3,000  rads  he  is 
placed  in  the  functionally  impaired  state  until  he  is  killed 
in  2  to  5  days. 

If  the  ACCOM. DOSE  is  greater  than  what,  he  can  tolerate 

(LETHAL. DCSE  (TROOP)  )  but  less  than  2,000  rads  the  TROOF  is 

unaffected  for  a  short  period  of  time  -  30  to  60  minutes. 

LET  T. IMPAIR. NOC  (TROCF)  =  R  ANDI.  F  (30 , 60 ,  1 ) 

At  the  end  cf  this  time  he  is  placed  in  the  functionally 

impaired  state.  He  is  killed  in  2  to  5  days. 

If  ACCOM. DOSE  is  less  than  LETHAL  .DOSE  (TP.OOP)  then  he  is 
fully  functional  and  will  not  be  killed  as  a  result  of  the 
current  detcnation. 
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K.  FUNCTION  0V2B.  PBESSUBE 


This  function  determines  the  distance  from  GZ  at  which  a 
target  element  will  encounter  a  50%  chance  of  survival  for 
the  achieved  7ield  and  required  oper pressure. 

I.  FUNCTION  DYNAMIC. EBESSO  BE 

This  function  determines  the  distance  from  GZ  at  which  a 
target  element  will  encounter  a  50 %  chance  of  survival  for 
the  achieved  yield  and  required  dynamic  pressure. 

B.  FUNCTION  ABBIVAL.TIME 

This  function  executes  the  arrival  time  function  and 
returns  the  time  after  detonation  when  the  blast  wave  will 
pass  the  given  distance  fromn  GZ. 

N.  FUNCTION  GAHMA 

Eguation  2.4  is  executed  by  this  function.  For  the 

achieved  yield  and  distance  from  GZ  the  amount  of  neutron 
radiation  is  returned. 

O.  FUNCTION  NEOTBON 

Eguation  2.3  is  executed  by  this  function.  For  the 

achieved  yield  and  distance  from  GZ  the  amount  of  neutron 
radiation  is  returned. 

P.  FUNCTION  BEAT 

Equation  2.2  is  executed  by  this  function.  For  the 

achieved  yield  and  thermal  requirement  the  distance  frcm  GZ 
where  a  50*  chance  of  lethal  burns  occurs  is  returned. 
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Q.  FOHCTICN  PRO BABI 1ITT 


This  function  approximates  the  area  under  a  N(0,1) 
probability  curve.  It  is  a  table  look-up  for  Fig.  4.8. 
Only  entry  arguments  (R/RD)  of  .5  and  greater  are  passed  to 
it.  If  the  argument  is  greater  than  1.5  the  probabilitiy  is 
set  tc  0. 


VI.  ENHANCEMENTS  AND  EXTENSIONS 


A.  OPTIMUM  ASSIGNMENT  OF  BATTEBIES. 

This  model  assigns  batteries  to  fire  on  an  essentially 
randoa  basis.  There  is  no  thought  given  to  optimizing  the 
coverage.  Choosing  the  highest  coverage  for  each  target  or. 
a  first  in  first  out  basis  may  not  provide  maximum  total 
coverage.  An  assignment  problem  must  be  solved  to  do  this 
correctly. 

Table  17  is  a  potential  nuclear  load  for  three  batteries 
in  the  following  example. 


T  ABLE  17 

Nuclear 

Load  by  Battery 

No.  of  rounds 

Battery 

Type 

.2  KT 

1  KT  2  KT  8 

1 

SRC 

1 

2  0  0 

2 

SRC 

1 

1  0  0 

3 

MRC 

0 

0  1  2 

Table  7  shoes  tfce  coverages  that  can  be  obtained  with 
the  various  battery/yieli  combinations  on  five  targets.  A 
zero  means  that  the  target  is  either  out  of  range  or 
coverage  is  less  than  command  guidance.  Note  that  battery  1 
has  2  rounds  of  1  KT  yield  and,  thus,  has  2  identical  rows 
for  the  1  KT  rounds.  A  similar  situation  is  true  for  the  8 
KT  rounds  of  battery  3. 
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TABLE  7 

Target  Coverage 


Target 

Battery  Yield  12345 


.2 

KT 

.72 

.00 

.00 

.00 

.00 

1 

1 

KT 

.82 

.80 

.00 

.45 

.00 

1 

KT 

.62 

.80 

.00 

.45 

.00 

.2 

KT 

.75 

.00 

.42 

.00 

.00 

2 

1 

KT 

.85 

.00 

.51 

.42 

.00 

2 

KT 

.92 

.63 

.47 

.  97 

.92 

< 

3 

KT 

.97 

.00 

.00 

.00 

.00 

8 

KT 

.57 

.00 

.00 

.00 

.00 

Tc  assign  battery/yield  combinations  to  a  target 
according  to  the  current  technique  would  produce  the 
following; 


Target  Eattery/Yield  Coverage 


1 

3/8 

Kt 

.97 

2 

1/1 

Kt 

.80 

3 

2/1 

Kt 

.51 

4 

3/2 

Kt 

.97 

This  arangement  results  in  a  total  coverage  cf 
target  number  5  cannot  be  engaged  because  the  only 
round  for  it  was  employed  on  target  number 
assignment  problem  that  needs  to  be  solved  is  one 
maximize  total  coverage.  This  case  is  very  easy  to 
inspection. 


3.25,  out 
feasible 
4.  The 
that  will 
solve  by 
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The  optimum  assignment  is; 

Target  Eattery/Yield  Coverage 

1  3/8  Kt  .97 

2  1/1  Kt  .80 

3  2/1  Kt  .51 

4  1/1  Kt  .45 

5  3/2  Kt  .92 

The  optimum  total  coverage  is  then  3.65. 

In  the  current  mcdel#  only  the  percentage  of  coverage  is 
used  to  elvaluata  the  value  of  the  target.  This  may  net  be 
the  item  to  maximize  if  the  various  targets  are  not  of  the 
same  type  and  population.  For  example,  a  97%  coverage  cn  an 
infantry  company  is  not  as  valuable  as  a  50%  coverage  on  a 
tank  regiment.  Tc  make  the  matrix  in  Table  v  mere 
realistic,  a  variable  such  as; 

value  =  coverage  +  element  value  *  target  population 
should  he  used. 


B.  VABI1TICN  OF  BSD 

Typically,  an  operations  order  specifies  the  trcop 
safety  guidance  as  a  negligible  risk  to  unwarned  exposed 
personnel.  This  is  the  category  used  in  the  medal  to 
determine  HSD.  Fig.  6.1  shows  the  different  separation 
distances  that  must  be  achieved  for  the  safety  categories 
cf ; 

1.  Unwarned  exposed. 

2.  Earned  exposed. 

3.  Earned  protected. 


Each  category  has  the  3  levels  of  radiation  exposure 
status  (RES)  of; 


1.  Negligible. 

2.  Moderate. 

3.  Emergency. 

UNWARNED  WARNED  WARNED 

EXPOSED  EXPOSED  PROTECTED 


RANGE 

NEG 

MOD 

EMER 

NEG 

MOD 

EMER 

NEG 

MOD 

EMER 

7500 

6900 

4500 

440  y 

4ino 

?2jr 

3400 

30  0  0 

1600 

JJOO 

7500 

0900 

4500 

4400 

«100 

33J0 

3500 

31  00 

1600 

4U  JO 

7500 

7000 

45jO 

44Jg 

4103 

.3300 

3S00 

3100 

1700 

■jOuO 

7t»C0 

70o0 

46J0 

4500 

4200 

3300 

36C0 

3700 

170  0 

oJjO 

7600 

7000 

4630 

4500 

4700 

.330*. 

3600 

3700 

17or. 

7u  jO 

/6C0 

7030 

46  JO 

4500 

4?00 

.3400 

3700 

3300 

looo 

cw  JO 

7600 

7100 

4600 

4500 

4200 

.34  J(. 

3700 

330  0 

lpno 

90  CO 

77C0 

71o0 

47  JO 

4bUU 

4  500 

340C 

3800 

3400 

looo 

i  uOoO 

7700 

71  JO 

47J0 

4600 

4  300 

3400 

39C0 

3500 

1POO 

ilO'jO 

7700 

7100 

47Jo 

4600 

4300 

.3500 

4000 

3*0  C 

looo 

:  tO'jQ 

7700 

7*00 

4700 

4600 

430(1 

3507 

4000 

3600 

1O00 

i  JujO 

7b00 

7iio0 

4800 

4700 

4400 

35  or 

4 1  CO 

3700 

1°00 

|4CU0 

7dtC 

7iu0 

4800 

4700 

4400 

3500 

4200 

37o0 

1°00 

iU0>,Q 

7bu0 

7*00 

4300 

4700 

4400 

36  Jl' 

4.300 

3or  o 

2000 

1  oUviO 

7n00 

7300 

4900 

4700 

4400 

36  J( 

4300 

3000 

7000 

»  7  JjQ 

7900 

7300 

49  Jo 

4800 

4900 

3600 

4400 

3°O0 

2P0C 

ldO.,0 

7^00 

73y0 

<9j0 

4800 

45Q0 

3hJ( 

4500 

4000 

2000 

Figure  6.1  Minim  urn  Separation  Distances. 


In  Fig.  6.1  the  MSD  for  a  negligible  risk  to  unwarned 
exposed  and  warned  exposed  personnel  are  7700  and  4600 
meters  respectively.  Recall  that  MSD  is  the  reason  for 
offset  and  offset  is  the  reason  for  some  very  severe 
degradation  in  coverage.  If  the  4600  meter  MSD  could  be 
used,  then  the  displacement  could  be  reduced  by  4300  meters. 
This  could  make  a  significant  difference  in  the  final 
outcome  of  the  simulation. 

Expansion  of  the  MSD  category  would  greatly  expand  the 
use  of  the  model  as  an  analytical  tool. 
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C.  COLLATERAL  DAMAGE  AVOID  AHCE 


Chapter  14  of  Ref.  2  contains  least  separation 
distances  (LSD)  to  avoid  damage  to  many  entities  such  as 
buildings,  forests,  bridges  and  civilian  populations.  This 
imposes  the  same  type  of  offset  requirement  on  the  aimpoint 
as  MSD  dees  except  that  it  precludes  damage  to  entities 
ether  than  friendly  troops. 

Use  of  this  idea  would  require  implementation  cf  a 


routine  similar  to  that  which  determines  offset  except  that 
the  description  of  what  is  being  avoided  must  also  be 
included. 

D.  HOITIELE  TARGETS 

As  it  has  already  been  shown,  the  center  of  a  target 
need  not  be  the  aimpeint.  It  would  be  highly  advantageous 
to  destroy  more  than  cne  target  with  a  single  burst.  Fig. 


6.2  shews  3  targets  and 


o 

Target  1 

Figure  6.  2 

displacement  distance  to  insure 
The  •'triangular'1  area  that  is 
feasible  area  for  each  target, 
anywhere  in  the  ''triangle”  will 


minimum  acceptable  coverage, 
inside  all  circles  is  a 
That  is,  a  round  employed 
destroy  all  3  targets. 
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A  method  to  evaluate  groups  of  targets  would  greatly 
expand  the  use  of  the  model  as  an  analytical  tool. 

E.  F ALLCUT 

Generally,  nuclear  weapons  are  detonated  at  or  above  an 
altitude  that  provides  a  99%  assurance  of  no  fallout  (Egs. 
3.4  and  3.5).  However,  fallout  may  be  an  effective  barrier 
to  movement  or  cause  mass  evacuation  of  a  contaminated  area. 
Control  cf  fallout  after  employment  is  impossible  due  to 
atmospheric  conditions  and  the  user  may  become  the  victim. 

Should  a  ground  burst  occur,  it  is  necessary  to  predict 
and  simulate  fallout.  Predicting  fallout  requires  that 
meterclogical  data  at  altitudes  up  to  and  including  cloud 
top  must  be  known.  A  deterministic  model  for  predicting 
fallout  is  relatively  simple  but  actually  simulating  fallout 
may  require  assistance  from  experts  in  meteorology. 


P.  ACTIONS  TAKEN  DOBING  AND  AFTEB  FALLODT. 

while  the  detonation  of  nuclear  weapons  will  certainly 
cause  seme  unique  actions  and  reactions,  there  are  some 
specific  actions  required  of  sddiers  upon  encountering 
fallcut. 

\m  Optimum  Time  of  Exit. 

For  personnel  caught  in  a  fallout  area  (presumably 
in  covered  fox  holes  cr  shelters)  there  is  an  optimum  time 
to  exit  the  area  to  minimize  the  absorption  of  nuclear 
radiation. 

While  soldiers  stay  in  the  protective  confines  of 
their  fox  holes,  they  take  advantage  of  the  shielding.  They 
also  prolong  their  exposure  time.  To  immediately  leave 
would  reduce  their  exposure  time  but  would  expose  them  to  a 
much  higher  dose  rate  outside  their  shelters. 
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For  given  initial  dose  rate,  decay  rate,  shielding 
and  exposure  time  to  a  safe  area,  there  is  an  optimum  time 
to  leave. 

2.  Crossing  Problem. 

Given  that  a  unit  must  cross  a  radiation  area,  there 
is  a  not-earlier-than  time  before  which  they  may  not  enter. 
It  is  a  function  of  EES,  fallout  pattern  and  decay  rate, 
speed  of  crossing  and  shielding  during  the  crossing. 
Inclusion  of  this  type  of  problem  into  the  the  maneuver 
section  of  a  driver  routine  would  certainly  have  an  impact 
upon  the  outcome  of  the  simulation. 

3 •  Decont  am inat ion. 

Actions  taken  to  decontaminate  vehicles  and 
personnel  are  relatively  simple  if  only  transportation  to 
and  from  a  decontamination  site  and  time  taken  to 
decontaminate  are  considered.  The  main  problem  with 
decontamination  is  one  of  supply. 


G.  ADDITIONAL  TARGET  ELEMENTS. 


The  current  model  has  only  tanks  and  troops  in  the 
target  areas.  Any  target  may  be  added  to  the  target  array 
by  defining  it  as  a  temporary  entity  and  allowing  a  target 
to  own  it.  To  damage  the  entity  requires  the  modeler  to 
call  one  or  more  of  the  damage  functions  with  the 
appropriate  arguments  and  then  dacide  if  the  entity  is 
damaged.  The  procedure  is  virtually  identical  to  what  is 
currently  used. 


APPENDIX  1 
BEGS ESS  ION  ANALYSIS 

A.  GANNA  B  ADI  AT  ION 

Since  the  acdel  needs  the  gamma  radiation  level  at  a 
specific  range  for  a  given  yield,  in  is  advantagous  to 
replot  Fig.  2.7  with  dose  as  the  dependent  variable.  Fig. 
A.  1  is  this  new  plot. 


Figure  A. 1  Gamma  Radiation. 

Surrendering  this  plot  to  the  APL  linear  regression 
package  fcr  the  equation; 

In  (Dose)  =  A  ♦  3  *  Range  ♦  C  *  In  (Yield) 

fails  tc  provide  an  adequate  solution.  Therefore,  the  model 
was  solved  in  pieces  and  reassembled. 


The  procedure  was  to  regress  each  line  independently  as 


In  <Dose)  =  A  ♦  B  *  Range 

The  coefficients  A  and  B  are  obtained  for 


for  each  yield.  The  coefficients  A 
each  yield  and  are  listed  in  Table  VI. 


TABLE 

VI 

Coefficients 

1  and  B 

Id 

A 

B 

11  .18 

-.0060 

11  .56 

-.0056 

11.76 

-.0054 

11.91 

-.0052 

12.02 

-.0052 

12.22 

-.0051 

12.22 

-.0050 

12.33 

-.0050 

12.38 

-.0049 

0 

12.42 

-.0049 

Fcr  each  yield,  the  coefficients  are  plotted  as  shown  in 
.  A. 2  and  A. 3. 


i. 


Figure  A.  3  Degression  Coefficient  B  (Data). 

The  coefficients  can  now  be  regressed  against  yield.  From 
Fig.  A. 2  the  first  ccefficient  appears  to  vary  as  the  square 
root  cf  yield.  Therefore,  the  equation; 

A  =  C0  ♦  Cj  •  /Yield 

is  fit  with  the  following  results. 

CB  *  10.74 


CT  =  .5599 


This  equation  is  plotted  in  Fig.  A. 4  along  with  the 
original  coefficient  data.  The  fit  is  quite  good. 


Figure  A. 4  Regression  coefficient  A  (Regressed) 
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There  is  a  temptation  to  say  that  the  coefficient,  E,  is 
really  a  constant  since  the  curve  is  very  flat  and  the 
represented  values  are  so  small.  Hovever,  it  must  be 

remembered  that  this  is  to  be  multiplied  by  the  range  and 
the  result  used  as  an  exponent.  Any  small  error  may  have 
catastrophic  results.  From  Fig.  A. 3  a  straight  line  of  the 
form; 

E  =  c0  ♦  C1  *  Yield 
is  fit  with  the  following  results. 

C0  =  -.0053 
C  =  .0001 

l 

The  equation  is  plotted  in  Fig.  A. 5  along  with  the  original 
coefficient  data. 


123456739  10 

Figure  A. 5  Regression  Coefficient  B  (Regressed). 

Assembling  the  equation  results  in  the  following 
equation. 


G  =  46166  • 


.5599/Yield  . 0001 -Yield *R 

e  e 

. 00 5  8  •  R 


(A.1) 


Equation  A.  1  is  plotted  on  Fig.  A. 6  for  even  yield 
along  with  the  original  data. 


l  1 
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Figure  A. 6  Full  Gamma  degression  Model. 

It  is  important  that  the  equation  be  accurate  in  the 
vicinity  cf  650  rads  as  this  is  the  mean  of  the  lethal  dose 
distribution  and  is  the  dividing  point  between  surviving  and 
dying. 


B.  BE0T5CN  BA  D I  AT  10  H 

Like  gamma  radiation,  the  model  needs  the  neutron 
radiation  level  at  a  specific  range  for  a  given  yield. 
However,  Fig.  2.6  appears  to  have  straight  lines  which  can 
te  exploited  in  a  regression  analysis.  Fig.  2.6  has  been 
converted  to  a  metric  scale  and  is  presented  in  Fig.  A. 7. 

The  neutron  regression  model  was  also  split  into  two 
submodels.  Each  line  was  regressed  independently  as; 


Range  =  A  ♦  B  *  In  (Yield) 


for  each  dose. 


The  coefficients  A  and  B  are  obtained  for  each  dose  and  are 
listed  in  Table  VII. 


TABLE  VII 
Coefficients 


Dose 

A 

B 

30 

1272 

171 

100 

1093 

158 

300 

9  18 

156 

1000 

766 

148 

30CC 

641 

129 

10000 

461 

138 

For  each  dose,  the  coefficients  are  plotted  as  shewn 
Figs.  A. 8  and  A.9.  The  coefficients  can  now  be  regres 
against  the  dose. 


For  the  first  coefficient.  A,  the  model  is 


In  (A)  =  C0  ♦  C  *  In  (Dose) 


30  100  300  1000  3000  10000 

Dose 

Figure  A.  9  Begression  Coefficient  B  (Data). 


For  the  second  coefficient,  B,  the  model  is 


with 


E  =  C0  ♦  C i  *  In  (Dose) 


C0  *  190.45 
Cj  =  -6.4145 

The  submodels  can  now  be  assembled  and  result  in  the 
following  full  model. 

Range  =  — -^yggg  +  (190-6 . 4  •  JLn(Dose  )  )  •  fcn(Yield)  (k  2) 
Dose’ 

What  the  model  needs  is  a  function  of  range  and  yield 
which  yields  dose.  Equation  A.  2  cannnot  be  separated  to  dc 
this,  therefore  an  approximation  is  offered.  The  term; 

190  -6.4  *  In  (Dose) 


is  the  slope  of  the  line  for  a  given  dose. 
30  to  10,000  rads,  the  average  slope  is  150 


For  doses  from 
The  proposal 


Range  (Meters) 


is  to  use  150  as  the  slope  for  every  line.  Fig.  A. 3  has  the 
data  plotted  again  as  well  as  Eg.  A. 2  and  the  new  proposal. 
For  doses  of  3,000  rads  and  below  the  new  proposal  seems  to 
be  an  iaprcveaent. 


Figure  A. 10  Heutron  Radiation  and  Regression. 


Table  VIH  lists  the  full  model  slope  and  the  percent 
error  if  a  slope  of  150  is  used.  Notice  that  the  errcr  i: 
negative  at  low  doses  but  is  getting  larger.  Thus,  it  i 
expected  to  be  in  perfect  agreement  at  some  dose. 


That  dcs 


C.  TBEBBAL  RADIATIOH 

The  model  uses  a  thermal  radiation  equation  which  yields 
range  as  a  function  of  yield  and  radiant  exposure.  Fig.  2.2 
is  already  in  that  format  with  the  exception  that  range  is 
in  miles.  Fig.  A.  1 1  is  the  same  figure  using  a  metric  range 
scale. 


3000 


rC*, 

c»>< 


«  1000 


a) 

to 

(2  500 


<  c>^  a 


Yield  (KT) 

Figure  A.  1 1  Thermal  Radiation. 

The  previously  used  method  of  using  submodels  will  again 
be  used.  The  full  model  is; 

Range  =  A  +  B  *  In  (Yield) 

for  each  exposure.  The  coefficients  A  and  3  are  obtained 

for  each  exposure  and  are  listed  in  tanle  IX  . 

For  each  exposure,  the  coefficients  are  plotted  as  shewn 
in  Figs.  A. 12  and  A. 13.  The  coefficients  can  now  be 
regressed  against  the  exposure. 
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TABLE  IX 

Regression  Coefficients  A  and  B 


Exposure  A  B 


3 

6.8331 

.5136 

5 

6.5  728 

.5051 

8 

6.3238 

.5239 

12 

6.1 137 

.5180 

25 

5.7  786 

.4926 

50 

5.4593 

.4773 

-r 


_ I _ I _ I _ I _ I - L_ 

3  5  3  12  2  5  50 

Cal/Cm2 

Figure  A.  12  Regression  Coefficient  A  (Data) 


For  the  first  coefficient.  A,  the  model  is 


A  *  C0  ♦  CT  *  In  (exposure) 


C  n  =  7.  3528 


Figure  A. 13  Regression  Coefficient  B  (Data). 


The  second  coefficient  could  now  be  regressed  against 
exposure,  but  several  models  tried  failed  to  yield  ar. 
accurate  explanation  of  behavior  of  the  coefficient  as  a 
function  of  exposure.  This  makes  sense  since  B  is  the  slope 
of  a  given  line  and  all  lines  appear  to  be  parallel. 
Therefore,  an  average  slope  value  is  \sed. 

B  =  .5051 

The  submodels  can  now  be  assembled  and  result  in  the 
following  full  acdel. 


Range  =  1561- (Yield) -5051  (i*4) 
(Exposure) -4885 

Normally,  a  figtre  would  be  used  to  compare  the  data  to 
the  regression  lines  at  this  point.  However,  the  agreement 
is  so  gocd  that  a  graph  would  lack  enough  detail  to  reveal 
any  differences.  Therefore,  Tables  X  and  XI  will  be  offered 
to  shew  hew  well  the  curve  fit  is. 


Original  Data 


pcsure 
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2 
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3 
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933 

1289 

1609 

2253 

2594 

2896 

c 

708 

1014 

1271 

1609 

1931 

2253 

8 

563 

805 

982 

1287 

1496 

1931 

12 
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644 
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1046 

1255 

1464 

25 

322 

451 

563 

724 

837 

998 

50 

24  1 

322 

386 
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TABLE  XI 

Hegression  Results 


Exposure 

1 

2 

Yield 

5 

7 

10 

3 

912 

1295 

1590 

2058 

2439 

2920 

5 

711 

1009 

1239 

1603 

1900 

2275 

8 

565 

802 

985 

1274 

1510 

1809 

12 

464 

658 

808 

1045 

1239 

1484 

25 

324 

460 

564 

730 

366 

1037 

50 

230 

3  28 

402 

521 

617 

739 

As  cne  can  clearly  see,  the  fit  is  excellent,  but  the 
saoe  caution  about  extrapolation  oust  apply.  Equation  A. 4 
is  valid  only  within  the  range  of  the  data. 


APPEHDII  B 
HOD EL  LISTIIG  FILE 


This  appendix  contains  a  cop;  of  the  listing  file  for  a 
test  run  of  the  ecdel.  There  are  three  targets,  five 
batteries  and  three  companies.  Each  target  has  50  tanks  and 
50  trcops  uniformly  distributed  within  the  ST.  Table  XII 
lists  the  target  data  while  Table  XIII  lists  the  battery 
data  and  Table  XIV  lists  the  company  data.  Target  types  are 
as  listed  cn  page  68. 


Target 

1 

Type 

TABLE  XII 

Target  Data 

X 

Y 

ST 

1 

Coordinate 

9500.0 

Coordinate 

0.0 

1400.0 

2 

4 

800  0.0 

0.0 

400.  0 

3 

5 

600  0.0 

0.0 

300.0 

TABLE  XIII 

Battery  Data 

Battery 

Size 

X 

T 

coordinate 

Coordinate 

1 

SRC 

0.0 

0.0 

i. 

HRC 

1000.0 

0.0 

3 

SRC 

1000.0 

500.0 

a 

HRC 

500.0 

0.0 

5 

SRC 

1500.0 

250.0 

TABLE  XIV 
Coapany  Data 


Coapany  X. 

coordinate 


1 

I 


200  0.0 
2000.0 
2000.0 


Coordinate 
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0.0 

-500.0 


Radius 
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APPENDIX  C 
GLOSSABY 

A.  BBAPCNS  EPPECTS 

Alpha  particle:  A  radioactive  emmision  identical  with  the 

helium  nucleus,  having  a  mass  of  four  and  an  electric  charge 
of  two  positive  units. 

Attenuation:  A  decrease  in  intensity  due  to  absorption  or 

scattering,  but  not  due  no  geometric  reduction  as  a  result 
of  range. 

Beta  particle:  A  radioactive  eoaision  identicle  in 

coopositicn  to  an  electron  with  the  possible  exception  that 
beta  particles  possess  a  very  high  speed. 

Induced  Badioactivity:  B adioactivity  produced  in  certain 

elements  by  the  capture  cf  free  neutrons.  In  the  area 
immediately  below  Ground  Zero  the  elements  of  sodium, 
■anganese,  aluminum  and  silicon  are  are  easily  induced. 
Neutrcn:  A  particle  of  mass  one  and  without  any  electric 

charge.  Neutrons  are  present  in  all  elements  except  light 

hydrcgen  and  are  required  to  initiate  the  detonation  process 
in  nuclear  weapons. 

Bad:  One  of  many  measures  of  radiation  density  or 

destructive  power.  It  represents  the  absorption  of  100  ergs 
of  ionizing  radiation  per  gram  of  absorbing  material,  such 
as  body  tissue. 

B.  TABGET  ANALYSIS 

Punctionally  impaired:  Personnel  who  exhibit  some  decreased 
ability  tc  perform  their  assigned  task,  but  are  not 
incapacitated.  Charicteriz  sd  by  vomiting,  diarrhea,  nausea, 
lethargy,  depression  and  mental  disorientation. 


T  .V.*  vr.  VJ  f 


T  •  - :  ■"■  ' ' 


Immediate  permanent  incapacitation  (IP)  Immediate  permanent 
incapacitation  is  the  ;most  severe  of  the  radiation  criteria 
and  causes  death  in  the  shortest  amount  of  time.  Personnel 
perforaing  tasks  which  are  physically  demanding  become 
incapacitated  within  5  minutes  of  exposure.  They  remain  in 
this  state  for  1  to  2  days  at  whick  time  death  occurs. 
Immediate  transient  incapacitation  (IT)  Personnel  receining 
this  dose  become  incapacitated  within  5  minutes  of  exposure 
and  remain  sc  for  30  to  45  minutes.  Personnel  then  recover, 
but  are  functionally  impaired  until  death  which  occurs  in  4 
to  6  days. 

latent  lethality  (LL)  Personnel  exposed  to  this  their  lethal 
dose  become  functionally  impaired  within  2  hours  of 
exposure.  Bore  than  half  of  this  group  will  die  within 
several  weeks. 

Incapacitated:  An  individual  who  performs  at  50  %  or  less 

of  his  pre-irradiation  performance  level. 

Badiation  exposure  status:  A  unit  wide  rating  relating  the 

average  accumulated  dose  to  the  risk  warranted  exposure  of 
pending  missions. 

Binimum  separation  distance  (BSD):  The  distance  friendly 

troops  must  be  from  the  burst  to  insure  less  than  2%  chance 
of  damage. 

Light  damage:  Light  damage  does  not  prevent  the  immediate 

use  of  an  itam.  Seme  repair  by  the  user  will  usually  be 
needed  to  make  the  item  operational. 

Hoderate  damage:  Hcderate  damage  prevents  use  of  an  item 

until  extensive  repairs  are  made.  This  degree  of  damage  is 
normally  sufficient  to  deny  use  of  the  equipment.  In  most 
situations,  achievement  of  this  degree  of  damage  will  be 
sufficient  to  support  tactical  operations. 

Severe  damage:  Severe  damage  prevents  use  of  an  item 

permanently.  There  may  be  some  situations  such  as  the 
attack  on  a  bridge  where  severe  damage  is  the  only  adequate 
degree  of  destruction. 
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